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Abstract 


Current procedures being used for the global prediction of 
300nm mean anomalies starting from known values of l°x 1° mean 
anomalies yield unreasonable prediction results when applied 
to 300nm blocks which have a rapidly varying gravity anomaly 
field and which contain relatively few observed-60nm blocks. 
Improvements to this situation are studied. 

Improvement of overall 300nm anomaly prediction is first 
achieved by using area-weighted as opposed to unweighted averaging 
of the 25 generated 60nm mean anomalies inside the 300nm block. 

Then, improvement of prediction over rough 300nm blocks is realized 
through the use of fully known l°x 1° mean elevations, taking 
advantage of the correlation that locally exists between 60nm 
mean anomalies and 60nm mean elevations inside the 300nm block. 

An improved prediction model which adapts itself to the roughness 
of the local anomaly field is found to be the model of Least 
Squares Collocation with systematic parameters, the systematic 
parameter being the slope b which is a type of Bouguer slope 
expressing the correlation that locally exists between 60nm mean 
anomalies and 60nm mean elevations. 
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1. Introduction 


Kaula (1966) described and implemented procedures for the 
prediction of mean gravity anomaly inside 300nm blocks (nm = 
nautical mile, see section 3). The procedures involved linear 
regression from known 60nm mean anomalies inside the individual 
300nm blocks. Rapp (1972) introduced two additions to the proce- 
dures and computer programs used by Kaula. The first, addition 
was the consideration of data noise by the introduction of a noise 
matrix into the prediction equations. The second addition was 
the estimation of accuracies of predicted 300nm anomalies. The 
techniques of Rapp actually involved the method of Least Squares 
Collocation with no systematic parameters. 

The procedures and computer programs used by Kaula, together 
with the subsequent enhancements by Rapp, constitute the current 
procedures being used to predict 300nm mean anomalies from l°x 1° 
mean anomalies. The purpose of this report is to study improvements 
to these current 300nm anomaly prediction procedures. 

First, background information which will be useful in subse- 
quent discussions will be covered in Sections 2 to 8. Then, 
in Section 9, it will be shown that current 300nm anomaly prediction 
procedures need to be improved over 300nm blocks which have a 
rough anomaly field and which contain relatively few observed 
60nm blocks. In Section 10, current procedures will be modified 
not necessarily to improve predictions over rough areas but to 
improve predictions through the use of more rigorous prediction 
formulas . 

To develop improvements to current prediction procedures 
as far as rough areas are concerned, three possibilities will 
be considered: 

1. the use of a local as opposed to a global covariance function 
in order that the local roughness of the anomaly field may be 
properly expressed in the prediction equations. This will be 
covered in Section 11. 

2. the removal of gross local trend between 60nm anomalies and 
elevations in accordance with the model 'Least Squares Collocation 
with systematic parameters.' This will be covered in Sections 

13, 14, 15, 16, and 17. 

3. the use of elevations as observed quantities, in addition 

to the observed anomalies, in a least squares collocation prediction 
involving an anomaly autocovariance function, an elevation auto- 
covariance function, and an anomaly-elevation crosscovariance 
function. This will be covered in Section 18. 
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2. Least Squares Collocation Equations 


Moritz (1972) gives a detailed presentation of Least Squares 
Collocation. Relevant equations are collected here for easy 
reference. 

Notations 


x 

A 

X 

AX 

s ' 
n 

E XX 
Cs ' s ' 
Cnn 
s 

^ss ' 

t 

B 

BX 

Css 

E tt 

E ss 


vector of observations 

a known rectangular matrix 

vector of systematic parameters 

systematic part of observations 

signal part of observations 

noise part of observations 

error covariance of estimated parameters 

covariance of signals of observation points 

covariance of noise of observations 

vector of signals at prediction points 

covariance between the signals being predicted and the 
signals at observation points 

vector of "complete signals" at computation points, defined 
through equation (6) 
a known rectangular matrix 

effect of the systematic parameters on the "complete signals" 
covariance of the signals at computation points 
error covariance of the predicted "complete signals" 
error covariance of the predicted signals 


Formulas 

The vector of observations is modeled as: 

x=AX+s'+n (1) 

The parameters are estimated from: 

X = (A t (C s . s , + Cnn) -1 A) -1 A T (C S , S , + C n n) _1 x . (2) 

The error covariance of estimated parameters is: 

E xx = <A T (C S . S , + Cnn)" 1 A)" 1 . (3) 

The vector of signals at computation points is predicted as: 

s = C gg i Q ( 4 ) 

where Q depends only on observation points and not on computation 
points. Q is given by: 

Q = (C s . s . + C nn ) _1 (x - AX). (5) 
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The vector of "complete signals" at computation points is defined as 


t = BX + s . (6) 

The error covariance of predicted "complete signals" is: 

= Cgg — C sg * (C g 1 g 1 + C nn ) 1 Cggi^ + (HA - B) ExxCA^hT — B^) (7) 

where 


H “ c ss'^ C s's' + C nn^ 1 • 

If it is assumed that the observation vector does not contain a 
systematic part, then the equations are: 


X = 

s 1 

' + n 

(9) 

t = 

s 

= ^ss'^s's 1 + C n n ^ x 

(10) 

E tt 

= 

E ss = c ss ~ Css-CCs's* + c nn) 1 c ss ,T 

(ID 


3. 300nm Block Subdivision Scheme 


The global 300nm block subdivision scheme employed by Kaula 
(1966) and retained in Rapp (1972) and in all predictions considered 
in this report is described in this section. The world is divided 
along parallels of latitude into 36 zones, each zone with a latitude 
extent of 5°. Individual latitude zones are then further separately 
divided along meridians such that near equal area blocks are 
obtained, with area as nearly equal as possible to the area of a 
5°x 5° block at the equator. The longitude extents of individual 
near equal area blocks are forced to be an integral number of 
degrees of longitude. Since an angular distance of 5° at the 
equator corresponds to a linear distance of 300 nautical miles, 
the near equal area blocks formed as described are called 300nm 
blocks . 

The specific algorithm used to generate the boundaries of 
300nm blocks is obtained after putting e equal to 5° in the 
equations found in Hajela (1975), pp. 2-3. Under this scheme 
longitude extents of 300nm blocks become 300nm± 30nm. Longitude 
extents of 300nm blocks vary from 5° of longitude at the equator 
to 120° of longitude at the poles. Figure 1 shows the appearance 
of boundaries of 300nm blocks in the Northern Hemisphere on an 
equal area projection. The boundaries in the Southern Hemisphere 
are symmetrically the same. For identification purposes the 300nm 
blocks of the world are uniquely numbered according to a wrap- 
around scheme, starting from the North Pole at 0°-longitude , 
proceeding eastward and southward as shown in Figure 1, until 
the last 300nm block at the South Pole is reached. The complete 
numbered set totalling 1654 blocks is shown in Figure 2. 
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E. 

A = 0° 


W. 



Figure 1 300nm Block Boundaries and 

Numbering Scheme 
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Each 300nm block is further conceptually divided into 25 
near equal area blocks. The 300nm block is first divided into 
5 l°-latitude zones. Each l°-latitude zone is then further divided 
along meridians ipto 5 near equal area blocks, called 60nm blocks, 
which have an integral number of degrees of longitude in longitude 
extent (see Hajela (1975), Table 4). Under this scheme the longi- 
tude extents of 60nm blocks are 60nm ±30nm. Longitude extents 
of 60nm blocks vary from one degree of longitude near the equator 
to 24 degrees of longitude near the poles. Within any 300nm block, 
the difference between two 60nm areas can reach, but never exceed, 
a factor of two (Hajela (1975) Table 4). 


4. Data Structure Inside the 300nm Block 


Procedures considered in this report predict a 300nm mean 
anomaly starting from known values of l°x 1° mean anomalies and 
l°x 1° mean elevations inside the 300nm block. Under the 300nm 
block subdivision scheme, there are considerably more l°x 1° blocks 
inside 300nm blocks near the poles than those near the equator. 

The number of l°x 1° blocks inside a 300nm blocks increases with 
increasing latitude of the 300nm blocks. For example, an equatorial 
300nm block has 5x5=25 l°x 1° blocks inside it. A polar 300nm 
block has 5x120=600 l°x 1° blocks inside it. Therefore, there 
are potentially much more l°x 1° data values inside 300nm blocks 
towards the poles than those towards the equator. 

To make the data distribution more uniform it is natural to 
average l°x 1° anomalies falling inside a particular 60nm block to 
obtain a single value to be considered the "observed" 60nm anomaly 
for that 60nm block. In the same way l°x 1° elevations inside the 
60nm block may be averaged to form the 60nm elevation. Since the 
component l°x 1° blocks of any 60nm block all lie in a common 1°- 
latitude zone, it follows that the areas of l°x 1° blocks inside 
a 60nm block are all equal; this makes it sufficient to straight 
average the l°x 1° values to form the 60nm value. The standard 
errors of the original l°x 1° values are propagated into the 60nm 
value according to the usual covariance propagation formula. The data 
structure is thus transformed such that the unit data blocks are 
near equal area 60nm blocks instead of l°x 1° blocks. The use of 
such a data structure is reasonable and greatly simplifies prediction 
computations . 

Therefore, during predictions the working data structure 
inside the 300nm block is as follows. There are exactly 25 60nm 
blocks arranged in a 5x5 array and numbered as shown in Figure 
3. There are 25 60nm mean anomalies, some or all of . which may 
be known, with each known 60nm anomaly having a standard error 
associated with it. There are 25 60nm mean elevations with cor- 
responding standard errors, all of which are always known for 
any 300nm block since l°x 1° elevations are fully known in the 
world. Finally, the exact areas of the individual 60nm blocks, 
useful for weighting purposes, are also known from the algorithm 
that generates 60nm block boundaries. 
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Anomaly predictions considered in this report will involve 
the predict ion of unknown 60nm mean anomalies inside the 300nm 
block (or, also possible, the direct prediction of the 300nm 
anomaly) using as known values, either just the known 60nm mean 
anomalies and their accuracies, or both the known 60nm anomalies 
and 60nm elevations and their accuracies. 
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Figure 3 


25 60nm blocks inside a 300nm block 


5. Covariance Matrices 


The computation of covariance matrices that enter the predic- 
tion equations considered in this report is described in this 
section. A numerical covariance function is first computed using 
global averaging procedures as described in Section 5.1. From 
this numerical covariance function the 25x25 covariance matrix 
of the 60nm blocks shown in Figure 3 can be formed as described 
in Section 5.2. This 25x25 covariance matrix becomes a basic 
covariance matrix from which the other covariance matrices needed 
in prediction equations are derived through covariance matrix 
propagation as discussed in Section 5.3. 


5.1 Computation of a 60nm Numerical Covariance Function 

Kaula (1966) and Rapp (1977) discuss a practical set of pro- 
cedures for the computation of an isotropic, homogeneous numerical 
covariance function for mean values. 

Let 


’ ^k : 

mean anomalies at two arbitrary blocks 

j and k 

> ^k • 

areas of blocks j and 

k 


^Jk : 

angular distance between 

the centers of 

blocks i and j 

*i : 

an angular distance. 
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Numerical mean anomaly covariances can be computed as follows: 


cov(Ag,Ag) = CCif^) = — A J A k Ag J Al l 

Z A j A k 


( 12 ) 


Z A. A k 


ip. = -k .yjk (13) 

1 ^ A. A k 

where the summation is taken over all blocks j and k satisfying: 

ijjjk = 0 for i=0 

0 < ih ^ < Aib for i=l (14) 



'"jk = 

0 

for 

o 

n 

•H 

0 < 


Aijj 

for 

i=l 

VI 

< 

ip .. < 

*jk 

i • Aip 

for 

i=2 , 3 


A\p is some angular distance that is considered approp- 
riate for use in the computations, considering the sizes of the 
blocks of interest. 

For 60nm blocks, Aib is chosen to be 

Aip = 1° (15) 

As explained in Rapp (1977), virtually identical covariances, 

out to ip approximately equal to 6°, are obtained even if 

the following simplifying condition is added to conditions (14): 

the 60nm blocks j and k fall inside a common 300nm block. 

(16) 

Using equations (12) and (13) along with conditions (14) 

(15) and (16), a 60nm numerical covariance function was computed, 
with the results shown in Table 1. The data set used was the 
OSU combined alt imeter /terrestrial data set of October 1979 (Rapp, 
1979, 1980). The anomalies of the data set were referred to the 
gravity formula of the Geodetic Reference System 1967. Also shown 
in Table 1 is the number of terms used in the summations occuring 
in equations (12) and (13). As can be seen, there is an ample amount 
of data for the determination of each C(^) . 

Given two 60nra blocks with angular distance between centers 
equal to i|i , the 60nm anomaly covariance between the two blocks 
can be interpolated from Table 1. In the procedures included 
in this report, linear interpolation is used; this corresponds 
to using a piecewise linear approximation to the covariance function 
(see Sunkel (1981a) for a discussion of implications of various 
interpolation methods). 
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Table 1 


cov(Ag,Ag) = C(^i>, computed from 
the combined altimeter /terrestrial data set 


i 


0 

1 

2 

3 

4 

5 

6 
7 


cov(Ag,Ag)=C(^j) No. of 
2 ^ Product 
Pairs Used 


(mgal 2 ) 


DEG 

1 RAD 

0.0000 

0.0000 

0.9796 

0.0171 

1.6201 

0.0283 

2.5081 

0.0438 

3.4661 

0.0605 

4.3769 

0.0764 

5.3623 

0.0936 

6.2280 

0.1087 


625 

37790 

382 

50057 

275 

92763 

207 

127132 

185 

101042 

176 

61813 

141 

8657 

109 

283 


5.2 Formation of Covariance Matrices 


From the discussed in Section 5.1 the 25x25 matrix 

containing the covariances among the 60nm blocks inside a 300nm 
block can be formed. This matrix will be denoted by C A A . 
Strictly, different 300nm blocks will have different CAgAg's 
simply because the set of ^-values for the component 60nm blocks 
will in general be different for different 300nm blocks. A possi- 
ble approach is therefore to take each 300nm block of the world, 
compute the ^-values referred to the centers of the component 
60nm blocks, and interpolate from the numerical covariance function 
for the covariances that will make up the C^g^g for the particular 
300nm block. However, this may take more computational effort 
than necessary. 

A simplifying assumption will instead be made that any 300nm 
block of the world is conceptually mapped into a common 300nm 
block shown in Figure 3; this common 300nm block has the properties 
that it is planar and that the angular distance ip between any 
two adjacent 60 nm blocks is considered to be equal to 1° along 
any row or column. Obviously, an equatorial 300nm block is distorted 
least in this kind of mapping while a polar 300nm block is distorted 
most. The distortions in positions translate into distortions 
of ^-values which in turn translate into distortions of the covar- 
iances interpolated from the numerical covariance function using 
the distorted ip -values. However, considering the well-known fact 
that the actual covariance values used do not critically affect 
the results of predictions (Rapp, 1977), coupled with the fact 
that 300nm blocks of greatest interest are hardly very near the 
poles, then the assumed mapping seems satisfactory. The advantage 
of such a simplifying assumption is in the great reduction of effort 
in the computation of covariances. 
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Referring to Figure 3, let 

cov(Ag^, Ag^): anomaly covariance between 60nm blocks i and j 

, Nj[ : row, column of block i 

Mj , Nj : row, column of block j . 

The distance between any two blocks i and j can now be computed 
in planar approximation as: 

\p = /(Mi-Mj) 2 + (Ni-Nj) 2 (degrees) . (17) 

The covariance cov(Agi , Agj) corresponding to this ^-value 
can be linearly interpolated from the numerical covariance function 
such as the one in Table 1. The 25x25 covariance matrix is then 
formed as: 

c AgAg = [cov(Ag ± , Agj ) ] (18) 

where i, j = 1,2,..., 25; and cov(Ag^,Agj) is the element 
of C AgA g in the ith row and jth column. 


5.3 Propagation of Covariance Matrices 

Let P and Q be two random vectors, by which will be 
meant two vectors which do not have any systematic part (M{P}=0, 
M{Q}=0). The covariance matrix between P and Q may be defined 
as : 

C PQ = M{PQ T } (19) 

where M is some averaging operator. 

Consider another random vector R which is related to 
P through the linear transformation: 


R 



( 20 ) 


where L Rp is the linear transformation matrix which transforms 
P to R . The covariance matrix between R and Q is: 


'RQ 


M {RQ t > 

M {L rp P Q t } 
= L rp M {P qt> 
= L Rp C P q . 


( 21 ) 


Again consider another random vector S which is related 
to Q through the linear transformation: 


s = l SQ Q (22) 

where LgQ is the linear transformation matrix which transforms 
Q to S . The covariance matrix between P and S is: 
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CPS = M {PS^} rp 
= M {PQT l| q } 

= M {PQT} t LJ q T 

C PS ” C PQ L SQ * 

(23) 

The covariance matrix between R and S is: 

C R o = M {RS t } 

= M { l rp p qT lsqTj 
= Lrp M {PQT) LsqT 
CrS = l RP c PQ LsQ 1 • 

(24) 

Equations (19), (20), (21), (22), (23), and 
how a covariance matrix CpQ propagates through 
formations of the random vectors P and Q . 

(24) define 
linear trans- 

Now let t and u be two non-random vectors, that is, 
vectors which contain both a systematic part and a random part 
(M{t}#0, M{u}^0 ) : 

t = t + T 

(25) 

u = u + U . 

(26) 

Here, an overbar denotes the systematic part, and a capital letter 
denotes the random part of a non-random vector. The covariance 
matrix between t and u is defined as: 

Cf u - M {(t-t) (u-u) T } 

(27) 

= M {TU T } 
C tu = C TU 

(28) 

Equation (28) says that the covariance matrix between two non- 
random vectors is equal to the covariance matrix between their 
random parts. 

Consider two non-random vectors v and w 


> 
+ 
i > 

ii 

> 


w = w + W 


Let v and w be related to each other through 
formation : 

the linear trans- 

v = Lw . 

(29) 

Then , 

v = L(w + W) 
= Lw + LW 


so that 

v = Lw 

(30) 

v ■ LW > -ii- 

(31) 


that is, if two non-random vectors are related to each other through 
some linear transformation matrix L , then their systematic and 
random parts are related through L . 

Using equations (28) and (31) it is now easy to extend covar- 
iance matrix propagation formulas to the case of non-random vectors. 
Let p , q , r , s be non-random vectors such that 


r = L rp p and s = L S q q . 
Then, by equation (31): 


R = Lrp P and S = L g Q Q . 

The covariance matrix between r and s becomes: 


C rs = Crs by equation (28) 

= Lrp CpQ LgQ T by propagation formula (24) 
= Lrp Cpq Lsq t by equation (28). 

Therefore , 


C 


rs 


L RP 


Cpq 



(32) 


Equation (32), analogous to equation (24), defines how a covariance 
matrix Cp q propagates through linear transformations of the 
non-random vectors p and q . 


From the basic covariance matrix c AgAg discussed in Section 
5.1, other covariances can now be derived through covariance matrix 
propagation. Let 


Ag 

Ag 


C ASAS 

C ASAS 


the 300nra anomaly 

vector with elements the 25 60nm anomalies inside the 
300nm block (Conceptually, there are always 25 60nm 
anomalies, though some of them may be unknown.) 
covariance between the 300nm anomaly and the 25 60nm 
anomalies 

variance of the 300nm anomaly. 


For the rest of this report a key assumption will be made 
as follows: 

If all the 25 60nm anomalies inside a 300nm block are known, 
then it is completely satisfactory to obtain the 300nm anomaly 
as the area-weighted average of those 25 60nm anomalies 
(60nm anomalies are obtained by straight averaging from l°x 1° 
anomalies as explained in Section 4). 

This assumption seems reasonable in practice; theoretical conse- 
quences may be pursued based on ideas contained in a recent report 
by Siinkel (1981b), but this will not be attempted at this time. 


- 12 - 


A direct consequence of the above assumption is that the 
300nm anomaly Ag will be related to the vector Ag containing 
the 25 conceptual 60nm anomalies through the linear transformation: 


Ag = VgAg ig - 


(33) 


where, specifically, L A -g Ag is the area-weighted averaging operator 
1 


L AgAg 


25 

I A ± 

i=l 


(Aj A 2 • • • A2 5 ) 


(34) 


with A, , the area of the 60nm block i , being known from the 
algorithm that generates 60nm block boundaries. Equation (34) is 
a generalization from current prediction procedures. Current 
prediction procedures assume that all the A-[ 's in equation (34) 
are equal, leading to a straight averaging in equation (33). 


Given equation (33) and applying covariance matrix propagation 
formula (21), 


C AgAg L AgAg °AgAg 


(35) 


That is, C A -g Ag is computable by area-weighted averaging elements 
of CAgAg along columns (or along rows, because of the symmetry 

of C AgAg } * 

Again considering equation (33) and applying propagation formula 
(23), T 

C A-gAg * L AgAg ■ <36) 

that is, the area-weighted average of the elements of C AgAg . 

Or, applying propagation formula (24), 

C A~gA"& = L A"gAg C AgAg L A'KAg (37) 

that is, the area-weighted average of all elements of c AgAg . 


To provide a feeling for the elements of C^Ag and c as as > 
the basic matrix CAgAg was first formed as described in Section 
5.1 and using numerical values from Table 1. Then, C Ag A g and 
c A"gAg were computed using equations (35) and (36) and assuming 
equal areas in equation (34). The results are shown in Table 
2. The values contained in C A g A g itself are not shown, these 
values being easily seen from equations (17) and (18) and Table 
1 • 
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Table 2 


Anomaly Covariance Between the 300nm Block and an 
Individual 60nm Block (i.e., cov(Ag , Agj ) ) , 
Assuming Equal Areas for the 60nm Blocks. 
Units: mgal 2 . 


226 

245 

250 

245 

226 

245 

267 

273 

267 

245 

250 

273 

279 

273 

250 

245 

267 

273 

267 

245 

226 

245 

250 

245 

226 


Variance of 300nm Block (i.e., cov(/Tg,ZTg) ) , Assuming Equal Areas 
for the 60nm Blocks: 252 mgal 2 . 


6. Equivalence of Indirect and Direct Predictions 


Let a vector t of "complete signals" be predicted from 
an observation vector x using the equations of Section 2. Now 
let another vector t be a vector of "complete signals": 

t = BX + s (38) 

where_ 

BX : systematic part of t 

s : signal part of t . 

Let t be related to t through the linear transformation: 

t = Lt (39) 


where L is the linear transformation matrix which transforms 
t to t . On substituting equation (6) into equation (39), 

t = L(BX + s) 

= LBX + Ls . (40) 

The systematic part of t becomes, on comparing equations (38) and 
(40) :_ 

BX = LBX . (41) 

The signal part of t becomes: 


s — Ls . 


(42) 
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Equations (38), (41), and (42) constitute an indirect prediction 
of t because in equation (42) the signal part of ^ is not 
predicted from observations directly but rather from a transformation 
of the predicted signal part of t . 

Substituting equation (4) into equation (42): 

i - L C ss , Q . (43) 

Considering the transformation equation (42) as being analogous 
to equation (20), and applying the covariance matrix propagation 
formula (21) : 

C- , = L C , . (44) 

ss' ss ' 

Therefore, equation (43) may be written as: 

i - C_ s , Q . (45) 

Equation (45) is analogous to equation (4). _Equations (38), (41), 
and (45) constitute a direct prediction of t because in equation 
(45) the signal part of t is predicted directly from observations. 


_The error covariance of the predicted "complete signals" 
in t is, using equation (7): 


E tt = C Ii - C S S ' <C s’s’ + C nn> _1 C ss ' 1 + 

+ (Ha - B) e xx (a t ht. b t > 


where 

H = C- , (C , , + C ) -1 

ss s 1 s ' nn 

From the covariance matrix propagation formula (23): 

T 

C— = L C L . 
ss ss 

Since equation (41) is an identity for all X , then 


(46) 

(47) 


(48) 


B = LB . 


(49) 


Substituting equations (44), (48), and (49) into 


Err = LC L - LC ,(C , . + C ) 

CS oc 1 0*0' nn ^ 


-1 


SS 


ss 


s ' s 


nn 


ss ' 


T l T 


+ L(HA - B)E xx (A T H T - B t ) L T 


ss ^ss'^s's' + ^nn^ 1 ■ + 


= L(C 
+ (HA - B) E 


ss 


XX 


(A T H t - B t ) ) l? 


(46) : 
+ 


Substituting equation (7) into the last equation, finally: 

E it = L E tt lT * < 50 > 

which also follows directly from equation (39) by covariance prop- 
agation . 
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Therefore, the_error covariance of the predicted "complete signals" 
contained in t may be computed using the direct method expressed 
by equation (46) or the indirect method expressed by equation 
(50) . 


The predictions included in this report use a combination 
of the indirect and direct methods described above. The indirect 
prediction is used to predict the 300nm anomaly; the t in equation 
(39) is the 300nm anomaly, and the t is the vector containing 
the 25 60nm anomalies inside the 300nm block. On the other hand, 
the direct method expressed by equation (46) is used to compute 
the variance of the predicted 300nm anomaly. 


7. Current Prediction Procedures 


Current prediction procedures to be examined in this report 
are those which predict a 300nm anomaly as the unweighted average 
of a completed set of 25 60nm anomalies inside the 300nm block 
(Kaula (1966), Rapp (1977)). The set of 25 60nm anomalies is 
completed by predicting any unknown 6 Onm anomalies from the known 
60nm anomalies inside the 30nm block using equations to be explained 
in Section 7.1. If there are no known 60nm anomalies inside the 
300nm block, the 300nm anomaly can still be predicted, using separate 
techniques which will not be examined in this report but which 
the reader can find in Rapp (1977). 

As part of current 300nm anomaly prediction procedures, a 
prior global analysis of the data set being used is done to generate 
a numerical covariance function. An example of a generated global 
numerical covariance function. is given in Table 1 (Section 5.1). 

The -numerical covariance function is used as explained in Section 
5.1 to form OigAg > the covariance matrix which is basic to 
the computation of the covariance matrices needed in the prediction 
equations . 


*7.1 Prediction of 60nm Anomaly 

To predict an unknown 60nm anomaly from known 60nm anomalies 
inside the 300nm block, Least Squares Collocation with no systematic 
parameters is employed by the current procedures. This means 
that an "observed" 60nm anomaly is simply modeled as being composed 
of a signal part and a noise part as in equation (9). The prediction 
equation when no systematic parameters are being modeled is equation 
(10), which with upgraded notations becomes: 

Ag i = C AgiAg*< C A£ftg* + D A|Ag^“ lA g* < 51 > 


where : 
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Agi 
A g* 
Ca giA g* 


C Ag* Ag* : 


predicted 60nm anomaly 

nxl vector of known 60nm anomalies 

covariance between the 60nm anomaly being predicted 
and the 60nm anomalies at observed 60nm blocks; lxn 
submatrix of the CAgAg discussed in Section 5.1 
covariance between the 60nm anomalies at observed 
60nm blocks; nxn submatrix of the C discussed 
in Section 5.1 

covariance of noise of known 60nm anomalies; diagonal 
matrix, with elements the square of the standard 
errors of the known 60nm anomaly values. 


7.2 Prediction of 300nm Anomaly 

Once the set of 25 60nm anomalies have been completed by 
predicting any unknown 60nm anomalies inside the 300nm block using 
equation (51), current procedures predict the 300nm anomaly as 
the straight average: 

25 

Ag = h> £ Ag i ' (52) 

i=l 

Comparing equation (52) with equations (33) and (34) it is seen 
that current procedures assume that all the 60nm areas A-^ are 
equal. The effect of this assumption will be examined in Section 
10 . 


Also, equation (52) is analogous to equation (42) so that 
current procedures use an indirect prediction of 300nm anomaly. 
However, for exact equivalence of equations (52) and (42) all 25 
Agj^ values in equation (52) and not just the unknown Agi 's must 
have been predicted by the methods of Section 7.1. Predicting 
the already known 60nm anomalies has the effect of filtering out 
the observational noise from those 60nm anomalies. The numerical 
effect of not filtering the known anomalies is mentioned in Section 
10 . 


The error variance of predicted Ag is computed by the direct 
method (Section 6) using equation (11), which with changed notations 
becomes : 

m Ag 2 = C AgAg “ C AgAg* (C Ag*Ag* + D AgAg ) 1 C AgAg* T (53) 

where 

C^-gAg* : covariance between the 300nm anomaly and the 60 nm 

anomalies at observed 60nm blocks; lxn subvector 
of the C^gAg given in equation (35) 

CAgAg : variance of the 300nm anomaly, given by equation (36) 

and ( 37) . 

Under the assumption of the current procedures that all the A^ 
in equation (34) are equal. 
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C A - gAg* : * xn subvec ' tor of c AgAg obtained by straight 
£j at5 averaging columns of C AgA g . See Table 2 for sample 

numerical values. 

C A - gA - : straight average of all elements of C A -g Ag or of 

C AgAg . See Table 2 for sample numerical value. 


8. The Test Area 


For numerical studies conducted in this report, use was made 
of 24 300nm test blocks covering the Tonga Trench area. The location 
of the 24 test blocks is shown in Figure 4. The numbers used to 
Identify the 300nm blocks follow the numbering scheme described in 
Section 3. The values shown inside the 300nm blocks arel°x 1° mean free 
air anomalies as contained in the OSU combined alt imeter/ terrestrial 
data set of October 1979 (the same data set used in Section 5.1). 

The axis of the Tonga lies along the highly negative anomaly values 
through blocks 1076, 1144, 1209, 1271, and 1272. 

For immediate comparison, Figure 5 shows l°x 1° mean free 
air anomalies for the test area as contained in the OSU terrestrial 
data set of October 1979 (Rapp, 1980). Where no value is indicated, 
there is no October 1979 terrestrial data for that particular 
l°x 1° block. As can be seen from Figures 4 and 5, the altimeter/ 
terrestrial data set is much more complete than the terrestrial 
data set. This is true not only in the test area but in the whole 
world, as can be seen from maps found in Rapp (1980) showing the 
global coverage of each of the two data sets. 


The test area is characterized by a wide range of roughness 
of 300nm anomaly field, from smooth to very rough. The specific 
character of each 300nm block regarding anomaly and elevation 
magnitude and roughness is given in Table 3. The definitions 
of descriptors used in Table 3 are given in Table 4, which uses 
the following abbreviations: 


ANOM 

RMSCA 

ELEV 

RMSCE 


300nm anomaly, computed by straight averaging known 60nm 
anomalies inside the 300nm block 

RMS centered anomaly, i.e., the root mean square of 60nm 
anomalies inside the 300nm block after the value ANOM 
has been subtracted from each 60nm anomaly 
300nm elevation, computed by straight averaging 60nm 
elevations inside the 300nm block. 

RMS centered elevation, i.e., the root mean square of 
60nm elevations inside the 300nm block after the value 
ELEV has been subtracted from each 60nm elevation. 
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Figure 4 : October 1979 l°x 1° Altimeter /Terrestrial 

Anomaly Coverage (mgals) and Location of 
300nm Test Blocks Over the Tonga Trench Area. 
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Figure 5 : October 1979 l°x 1° Terrestrial Anomaly Coverage 

(mgals) and Location of 300nm Test Blocks 
Over the Tonga Trench Area 
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Table 3 


Characterization of the 24 300nm Test Blocks 
Covering the Tonga Trench Area, Regarding 
Anomaly and Elevation Magnitude and Roughness 


Anomaly Elevation 


Block 

Magnitude 

Roughness 

Magnitude 

Roughness 

1005 

MED HIGH 

MILD 

MED LOW 

SMOOTH 

1006 

MEDIUM 

MILD 

LOW 

MED RGH 

1007 

MEDIUM 

MILD 

LOW 

MED RGH 

1008 

MEDIUM 

MILD 

LOW 

MED RGH 

1075 

MED HIGH 

SMOOTH 

MED LOW 

MED RGH 

1076* 

MED HIGH 

VRY RGH 

LOW 

VRY RGH 

1077 

MEDIUM 

MILD 

VRY LOW 

SMOOTH 

1078 

MEDIUM 

SMOOTH 

LOW 

SMOOTH 

1142 

MED HIGH 

MILD 

LOW 

MED RGH 

1143 

MED HIGH 

MILD 

MED LOW 

MED RGH 

1144* 

MEDIUM 

VRY RGH 

LOW 

VRY RGH 

1145 

MEDIUM 

SMOOTH 

VRY LOW 

SMOOTH 

1208 

MED HIGH 

SMOOTH 

LOW 

MED RGH 

1209* 

MED HIGH 

VRY RGH 

LOW 

VRY RGH 

1210 

MED HIGH 

MILD 

VRY LOW 

SMOOTH 

1211 

MEDIUM 

SMOOTH 

VRY LOW 

SMOOTH 

1270 

MED HIGH 

SMOOTH 

MED LOW 

RGH 

1271* 

MEDIUM 

VRY RGH 

LOW 

RGH 

1272* 

MED HIGH 

RGH 

VRY LOW 

SMOOTH 

1273 

MEDIUM 

SMOOTH 

VRY LOW 

SMOOTH 

1329 

MED HIGH 

MILD 

MED LOW 

MED RGH 

1330* 

MEDIUM 

RGH 

LOW 

RGH 

1331 

MEDIUM 

SMOOTH 

VRY LOW 

SMOOTH 

1332 

MEDIUM 

SMOOTH 

VRY LOW 

SMOOTH 

♦rough to 
in this 

very rough 
report . 

300nm blocks, 

critical to the 

studies 
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Table 4 


Definitions of Descriptors for 300nm Blocks 


A. Anomaly Magnitude. ANOM = Anomaly 

% Out of 

Descriptor Range of Values Totall 


LOW ANOM < -30 mgals 4 
MED LOW -30 < ANOM < -10 20 
MEDIUM -10 < ANOM < 10 51 
MED HIGH 10 < ANOM < 30 21 
HIGH 30 < ANOM 4 


B. Anomaly Roughness. RMSCA = RMS Centered Anomaly. 


SMOOTH 0 < RMSCA < 10 mgals 42 
MILD 10 < RMSCA < 25 43 
RGH* 25 < RMSCA < 45 11 
VRY RGH* 45 < RMSCA 4 

C. Elevation Magnitude. ELEV = Elevation. 

VRY LOW ELEV < -5200 m 6 
LOW -5200 £ ELEV < -3000 46 
MED LOW -3000 5 ELEV <0 19 
MEDIUM 0 1 ELEV < 500 17 
HIGH 500 < ELEV < 1200 8 
VRY HIGH 1200 < ELEV 4 


D. Elevation Roughness. RMSCA = RMS Centered Elevation. 


SMOOTH RMSCA < 500 m 62 
MED RGH 500 < RMSCA < 1000 20 
RGH 1000 < RMSCA < 1700 13 
VRY RGH 1700 < RMSCA 5 


♦critical categories for the purpose of this report. 
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Table 4 was established from an analysis of the combined 
altimeter/terrestrial data set. Since this data set had adequate 
l°x 1° mean anomaly coverage (Rapp, 1980) and since l°x 1° mean 
elevations are fully known, the computations of the quantities 
ANOM, RMSCA , ELEV, and RMSCE for the 300nm blocks were generally 
meaningful. Each 300nm block of the world was taken, and values 
of ANOM, RMSCA, ELEV, and RMSCE computed for the block. These 
values then fell within corresponding ranges given in Table 
4. A frequency count was kept of the number of 300nm blocks 
falling within the various ranges, and the percentages out 
of the total number of blocks considered were computed and 
given in the last column of Table 4. For example, Table 4 
says that 300nm blocks which have RMSCA ^ 45 constitute only 
4% of the total number of observed 300nm blocks; as such, 
blocks with RMSCA ;> 45 may be regarded as having very rough 
anomaly fields. Similary, 300nm blocks with 25 1 RMSCA < 45 are 
seen to constitute about 11 l of the total number of 300nm blocks 
considered, and these may be regarded as blocks with rough 
anomaly fields. It will be seen in Section 9 that 300nm blocks 
with rough to very rough anomaly fields are the blocks where 
current 300nm anomaly prediction procedures need to be improved. 

In Table 3, the blocks marked with asterisks are critical 
blocks for studies conducted in this report (see Section 9). 

These blocks fall in the categories of blocks with rough to 
very rough anomaly fields, the categories marked also with 
asterisks in Table 4. It should be noted from Table 3 that 
the blocks categorized as having rough to very rough anomaly 
fields are the same blocks which are categorized as having 
rough to very rough terrain, a consequence of anomaly-elevation 
correlation that will be more fully covered in Section 12. 


9. Problem With Current Procedures 


To examine the performance of current prediction procedures 
under a variety of cases of anomaly field roughness and number 
of known 60nm blocks, the procedures were used to predict 300nm 
anomalies for the 24 test blocks in the Tonga Trench area. 

Since all the test blocks have their 25 60nm anomalies known 
in the altimeter /terrestrial data set, "true" values of the 
300nm anomalies were obtained by area-weighted averaging the 
25 known 60nm anomalies inside the individual 300nm blocks, 
in accordance with equation (33). 

Then, to simulate lack of data during prediction, only 
those l°x 1° altimeter /terrestrial anomalies for which there 
was a corresponding known l°x 1° anomaly in the terrestrial 
data set were used to predict 300nm anomalies. Therefore, 
the data configuration was that shown in Figure 5, although 
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the actual data values used came from the altimeter /terrestrial 
data set shown in Figure 4. Predictions were done using current 
procedures (Section 7) and the global numerical covariance 
function given in Table 1. The results are shown in Table 
5. 


Column 1 of Table 5 gives the standard block sequence 
number. Column 2 gives N , the number of known 60nm anomalies 
used in the prediction. Column 2 also gives the character 
of the anomaly field as to roughness, obtained from Table 3. 
Column 3 gives "true" values. Column 4 gives the predicted 
300nm anomaly. Column 5 gives the error of prediction (predicted 
value minus "true" value) and the standard error of prediction. 

The immediate conclusion from Table 5 is that current 
prediction procedures yield satisfactory results for smooth 
to mild 300nm anomaly fields, start to have problems with 
rough 300nm anomaly fields, and definitely yield unreasonable 
results for very rough 300nm anomaly fields. As Table 4 
shows, rough 300nm blocks (25 < RMSCA < 45) constitute about 
10% and very rough blocks (RMSCA k 45) about 5% of all the 
300nm blocks of the world, so that roughly 15% of all 300nm 
blocks are problematic blocks similar to blocks 1076, 1144, 

1209, 1271, 1272, and 1330 shown in Table 5. Other tests 
would indicate that the problem is particularly pronounced 
for rough to very rough 300nm blocks which contain relatively 
few known 60nm blocks, say less than 15, and is rather insig- 
nificant when the number of known 60nm blocks is close to 
the maximum number 25, say greater than 20. 


10. Prediction With Area-Weighted Averaging 


A 300nm anomaly prediction procedure which is more rigorous 
than the current procedrues may immediately be obtained by 
using area-weighted averaging expressed by equations (33) and 
(34) instead of the straight averaging expressed by equation 
(52). Also, pursuing the second comment under equation (52) 
all 25 60nm anomalies inside the 300nm block, unknown or known, 
must be predicted prior to area-weighted averaging the 25 60nm 
anomalies. 

Therefore, prediction of the 300nra test block anomalies 
was done with what will be called Method A. The data set used 
was the October 1979 terrestrial data set. The covariance 
function used was that given in Table 1, which was a covariance 
function based on the altimeter /terrestrial data set. Strictly, 
a new covariance function should have been computed based on 
the terrestrial data set, but this would not have made any 
difference for the purpose of this section. 
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Table 5 


Prediction Using Current Procedures With 
Global Anomaly Covariance Function, 
Illustrating Prediction Problems 
Over Rough 300nm Blocks 


Block 

N/CHAR 

Anomaly Error/m^g 

(mgal) (mgal) 

"True" Current Proc Current Proc 
Anom Alt /Ter Alt/Ter 

(mgal) Subset Subset 

1005 

16/MLD 

11.1 

8.8 

-2. 3/3.5 

1006 

16/MLD 

9.2 

6.3 

-2. 9/3. 2 

1007 

22/MLD 

2.4 

0.5 

-1. 9/2.0 

1008 

16/MLD 

-2.9 

-2.1 

0. 8/3.1 

1075 

23/SMTH 

19.4 

20.4 

1.0/1. 8 

1076* 

13/VRY RGH 

21.7 

-3.7 

-25.4/5.6 

1077 

20/MLD 

5.3 

7.1 

1.8/2. 5 

1078 

13/SMTH 

-4.3 

-3.9 

0. 4/3.9 

1142 

11/MLD 

23.4 

21.4 

-2. 0/5. 5 

1143 

3/MLD 

21.6 

13.8 

-7.8/11.0 

1144* 

22/VRY RGH 

10.3 

5.5 

-4. 8/2.0 

1145 

8/SMTH 

4.1 

6.0 

1. 9/6.1 

1208 

3/SMTH 

16.5 

8.4 

-8.1/10.8 

1209* 

13/VRY RGH 

16.0 

-8.8 

-24.8/3.6 

1210 

10/MLD 

17.9 

10.0 

-7. 9/6. 3 

1211 

ft 

-0.5 

ft 

ft 

1270 

2/SMTH 

22.2 

19.0 

-3.2/11.8 

1271* 

13/VRY RGH 

17.0 

-3.6 

-20.6/5.9 

1272* 

5 /RGH 

13.4 

-4.5 

-17.9/9.4 

1273 

7/SMTH 

0.7 

-0.2 

-0.9/7. 7 

1329 

19/MLD 

27.1 

26 . 5 

-0. 6/3.0 

1330* 

21 /RGH 

2.4 

-0.7 

-3. 1/2. 4 

1331 

ft 

7.7 

ft 

ft 

1332 

8/SMTH 

6.2 

0.0 

-6. 2/7. 9 


#No October 1979 terrestrial data exists. 

♦blocks showing prediction problems 
(rough and very rough blocks) 
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Method A used previous equations summarized as follows. 
First, the 25 60nm anomalies were predicted using an extension 
of equation (51): 


Ag = C 


AgAg* 


(C 


Ag*Ag* 


D a * a 
A gAg 


-l 


A g* 


(54) 


where now Ag is a 25 x 1 vector containing the 25 predicted 
60nm anomalies. Then, using the area-weighted averaging operator 
given by equation (34): 


L AgAg 


1 

25 

_.L A i 


(Ai A 2 


i=l 


A 2 5 ) 


the 300nm anomaly was indirectly predicted (Section 6) using 
equation (33): 


A~g = 


L AgAg 


Ag 


Finally, the variance of the predicted 300nm anomaly was computed 
using equation (53): 

m A~g 2 = C A"gAg " C A~gAg* (C Ag*Ag* + D A|Ag° * C AgAg* 


where C£% A g* and c A "g A ~g were computed by area-weighted averaging 
from the basic covariance matrix C A g Ag . The prediction results 
and the square root of m^-g 2 are shown in Table 6, column 3. 

A direct prediction procedure which is exactly equivalent 
to the indirect procedure used by Method A may be done as 
follows, if desired. Equation (54) when substituted into equation 
(33) results in: 


Ag ” L AgAg C A gAg* (C 


Ag*Ag* 


+ D. * . 
AgAg 


' Ag* 


Applying covariance matrix propagation, the last equation may 
be written in the form of the direct prediction: 

Ag = C Agi g* (C AgV + D Ag*4g* rl ’ (55> 

The variance of predicted Ag is again given by equation (53). 

For the evaluation of Method A, current prediction procedures 
were applied to the same terrestrial data set used by Method 
A. The covariance function used was again the one given in 
Table 1. The results of predictions are shown in Table 6, 
column 2. The difference between Method A and the current 
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Table 6 


Predictions Using Method A With Global Anomaly 
Covariance Function, Illustrating the Effect of 
Area-Weighted Averaging and Prediction of 
all 25 60nm Anomalies. 


Block 

Anom/m^g 
(mgal ) 

Current Proc 
0ct79 Terr Set 

O 

Anom/m^g 
(mgal) 
Method A 
0ct79 Terr Set 

Anom 

Difference 
(mgal) 
Method A - 
Current Proc 

1005 

3.4/5. 5 

4. 8/5. 5 

1.4 

1006 

13.0/5.3 

12.3/5.3 

-0.7 

1007 

7. 9/3. 8 

8. 3/3. 8 

0.4 

1008 

8.9/5. 3 

5. 5/5. 3 

-3.4 

1075 

21.8/3.8 

18.6/3.7 

-3.2 

1076* 

3. 7/7. 5 

11.8/7.5 

8.1 

1077 

2. 4/4. 5 

3. 7/4. 5 

1.3 

1078 

-4. 9/5. 6 

-4. 6/5. 6 

0.3 

1142 

34.8/7.4 

32.4/7.4 

-2.4 

1143 

14.1/12.3 

13.5/12.4 

-0.6 

1144* 

-1.3/4. 5 

1.6/4. 5 

2.9 

1145 

-15.9/8.5 

-14.1/8.5 

1.8 

1208 

0.6/12.4 

0.5/12.8 

-0.1 

1209* 

-20.3/7.6 

-6.6/ 7.8 

13.7 

1210 

-11.1/8.0 

-10.1/8.0 

1.0 

1211 

ft 

If 

If 

1270 

25.7/13.2 

23.5/13.4 

-2.2 

1271* 

-3. 2/7. 6 

0.4/7. 5 

3.6 

1272* 

-17.1/11.1 

-13.5/11.4 

3.6 

1273 

-16.0/8.7 

-15.4/9.0 

0.6 

1329 

15.3/5.4 

27.4/5.1 

12.1 

1330* 

-13.2/5.0 

-15.1/4.8 

-1.9 

1331 

If 

ft 

If 

1332 

-23.8/10.5 

-20.9/10.7 

2.9 

# No October 

1979 terrestrial 

data exists. 



♦rough and very rough blocks. 
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procedures is indicated in Table 6, column 4. It is seen from 
Table 6 that there tend to be cases of considerable difference- 1 - 
between the two methods in their predicted 300nm anomalies, 
although the two methods yield practically the same predicted 
accuracies. Since Method A is more rigorous than current proce- 
dures, Method A has to be the preferable method. 

To study the numerical effect of not predicting the already 
known 60nm anomalies, area-weighted averaging (Method A) was 
again applied, but this time without predicting the known 60nm 
anomalies inside the 300nm block. The same data set and covariance 
function as used before were utilized. The method is called 
Method A' and the results are shown in Table 7, column 3. The 
difference between Method A' and Method A is shown in Table 

7, column 4. It is seen that there tend to be cases of consid- 
erable difference between the two methods. It was decided that 
it would be worth the effort to predict all 25 60nm anomalies 
in an improved method for the prediction of 300nm anomaly. 


As a further test of Method A, the method was applied 
to a subset of the altimeter/terrestrial data set obtained 
by considering only those l°x 1° altimeter /terrestrial anomalies 
for which there were corresponding known l°x 1° anomalies in 
the terrestrial data set. The covariance function used in 
the predictions was again the one given in Table 1. The results 
which are given in Table 8 again show that there tend to be 
cases of considerable difference between Method A and current 
procedures, this time in the case of altimeter/terrestrial 
data. An important case occurs for block 1271 for which Method 
A gave a much better 300nm anomaly prediction than current 
procedures judging from closeness of the predicted value to 
"true" value given in column 2 of Table 8. The smaller values 
of predicted accuracies in Table 8, column 4 as compared with 
the predicted accuracies in Table 6, column 3 are due to the 
smaller standard deviations in the altimeter/terrestrial data 
as compared with those in the terrestrial data. 


This report will often deal with comparison of alternative 
prediction methods in terms of differences in predicted results. 
Whether a certain difference is large or small should always 
be decided in the light of the magnitude of the estimated predic- 
tion error. A difference of lOmgals may be considered very 
large when the estimated prediction error over the 300nm block 
is just 3mgals, but it is certainly not that large when the 
estimated prediction error is 15mgals. Regarding the magnitude 
of the estimated prediction error, an allowance should be made 
for the fact that the use of a global covariance function in 
the prediction will tend to give a too small value for the 
estimated prediction error over a rough area. 
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Table 7 


Predictions Using Method A' With Global Anomaly 
Covariance Function, Illustrating the Effect of 
Not Predicting all 25 60nm Anomalies in Method A 



Anom/m^-g 

Anomaly 

Anom 


(mgal ) 

(mgal) 

Difference 


Method A 

Method A' 

(mgal) 


Terrestrial 

Terrestrial 

Method A'- 

Block 

Data Set 

Data Set 

Method A 

1005 

4. 8/5. 5 

3.3 

-1.5 

1006 

12.3/5.3 

13.0 

0.7 

1007 

8. 3/3. 8 

7.9 

-0.4 

1008 

5.5/5. 3 

8.9 

3.4 

1075 

18.6/3.7 

20.7 

2.1 

1076* 

11.8/7.5 

3.7 

-8.1 

1077 

3. 7/4.5 

2.3 

-1.4 

1078 

-4. 6/5. 6 

-4.9 

-0.3 

1142 

32.4/7.4 

34.9 

2.5 

1143 

13.5/12.4 

14.5 

1.0 

1144* 

1.6/4. 5 

-1.1 

-2.7 

1145 

-14.1/8.5 

-15.9 

-1.8 

1208 

0.5/12.8 

0.5 

0.0 

1209* 

-6. 6/7. 8 

-23.3 

-16.7 

1210 

-10.1/8.0 

-11.1 

-1.0 

1211 

# 

ft 

ft 

1270 

23.5/13.4 

27.1 

3.6 

1271* 

0.4/7. 5 

7.0 

6 ,6 

1272* 

-13.5/11.4 

-15.8 

-2.3 

1273 

-15.4/9.0 

-15.9 

-0.5 

1329 

27.4/5.1 

17.4 

-10.0 

1330* 

-15.1/4.8 

-15.0 

0.1 

1331 

# 

ft 

# 

1332 

-20.9/10.7 

00 

• 

CO 

<N 

1 

-2.9 

#No October 

1979 terrestrial 

data exists. 



♦rough and very rough blocks. 
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Table 8 


Predictions Using Method A With Global Anomaly 
Covariance Function, Applied to the Combined 
Altimeter /Terrestrial Subset 


Block 

"True" 
Anom 
(mgal ) 
Alt /Ter 
Subset 

Anom/m^g 
(mgal) 
Current 
Procedures 
Alt/Ter Subset 

Anom/mA'cr 
(mgal) 
Method A 
Alt /Ter Subset 

Anomaly Diff 
(mgal) 
Method A- 
Current 
Procedures 

1005 

11.1 

8. 8/3. 5 

8. 7/3. 5 

-0.1 

1006 

9.2 

6.3/3. 2 

6. 5/3. 2 

0.2 

1007 

2.4 

0. 5/2.0 

0. 1/2.0 

-0.4 

1008 

-2.9 

-2. 1/3.1 

-2. 3/3.1 

-0.2 

1075 

19.4 

20.4/1.8 

19.0/1.6 

-1.4 

1076* 

21.7 

-3. 7/5. 6 

-4. 1/5. 5 

-0.4 

1077 

5.3 

7.1/2. 5 

7.0/2. 5 

-0.1 

1078 

-4.3 

-3. 9/3. 9 

-4. 0/3. 9 

-0.1 

1142 

23.4 

21.4/5.5 

21.3/5.4 

-0.1 

1143 

21.6 

13.8/11.0 

13.3/11.0 

-0.5 

1144* 

10.3 

5. 5/2.0 

6. 0/2.0 

0.5 

1145 

4.1 

6. 0/6.1 

6. 0/6.1 

0.0 

1208 

16.5 

8.4/10.8 

8.2/11.3 

-0.2 

1209* 

16.0 

-8. 8/3. 6 

-11.9/3.7 

-3.1 

1210 

17.9 

10.0/6.3 

10.0/6.3 

0.0 

1211 

-0.5 

tt 

ft 

H 

1270 

22.2 

19.0/11.8 

19.0/11.8 

0.0 

1271* 

17.0 

-3. 6/5. 9 

7. 3/5. 7 

10.9 

1272* 

13.4 

-4. 5/9. 4 

-2. 7/9. 7 

1.8 

1273 

0.7 

-0.2/7. 7 

-0. 2/8.0 

0.0 

1329 

27.1 

26.5/3.0 

27.9/2.7 

1.4 

1330* 

2.4 

-0.7/2. 4 

-4. 3/2.1 

-3.6 

1331 

7.7 

H 

# 

# 

1332 

6.2 

0.0/7. 9 

-0.2/8. 1 

-0.2 


#No October 1979 terrestrial data exists, 
♦rough and very rough blocks. 
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In summary, Tables 6, 7, and 8 indicate that the two ideas 
behind Method A are worth incorporating in an improved 300nm 
anomaly prediction. The two ideas are: 

1. prediction of all the 25 60nm anomalies, unknown or known, 
inside the 300nm block. This has the effect of filtering 
out observational noise from the known 60nm anomalies; 

2. Use of area-weighted averaging of the 25 60nm anomalies 

to obtain the 300nm anomaly. This will improve prediction 
results especially for rough 300nm blocks with greatly 
unequal 60nm areas inside them (as mentioned in Section 
3, the difference between two 60nm areas inside a 300nm 
block can reach, though cannot exceed, a factor of two). 
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11. Prediction With Rough Covariance Function 


Covariance matrices used in Least Squares Collocation equations 
are formed from a basic covariance function as described in Section 
5.2. A covariance function such as the one given in Table 1 is 
a global covariance function, having been obtained by global averaging 
from 300nm blocks both rough and smooth. A global covariance 
function cannot truly represent the local anomaly behavior over 
rough 300nm blocks because rough 300nm blocks form a small minority 
(see Table 4) out of the total number of blocks used in the averaging 
process to obtain the global covariance function. It is therefore 
interesting to investigate whether 300nm anomaly predictions over 
rough areas can be improved through the use of a rough covariance 
function obtained by averaging from only those 300nm blocks with 
rough anomaly field. 

Table 9 shows a rough numerical covariance function cov(Ag,Ag). 
The function was computed from the combined altimeter/terrestrial 
data set using the same procedures ued to generate Table 1. This 
time, only those 300nm blocks which were classified as very rough. 

(see Table 4, RMSCA > 45mgals) were used in the averaging process 
that computed cov(Ag,Ag). Note the big difference between the 
covariances from Tables 1 and 9; in particular, note the very 
much enlarged variance (the covariance at ip=0°) and the shortened 
correlation length (the ip -value at which the covariance falls 
to half the value of the variance). 


Table 9 

cov(Ag,Ag) from the OSU Combined Altimeter /Terrestrial Data Set 
of October 1979, Using only Those 300nm Blocks 
With RMSCA > 45mgals (see Table 4) 


cov( Ag, Ag) No. of Product 

Deg. Rad. (mgal 2 ) Pairs Used 


0.0000 

0.0000 

4222 




1470 

0.9793 

0.0171 

1934 




1986 

1.6247 

0.0284 

769 




3620 

2.5074 

0.0438 

84 




5029 

3.4714 

0.0606 

17 




4093 

4.3831 

0.0765 

66 




2495 

5.3817 

0.0939 

-229 




353 

6.2559 

0.1092 

-1797 




10 

The rough covariance 

function 

from Table 9 

was used as described 

in Section 5.2 to form covariance matrices 

needed in collocation 

equations. 

Subsequently, 

predictions 

were 

done 

under current 

procedures as described in 

Section 

7, 

this 

time 

using the rough 

covariance matrices just formed. Prediction results are given 

in Table 10 , 

Column 3. For 

comparison , 

results of predictions 


- 32 - 


Table 10 


Predictions Illustrating the Use of a Rough Covariance Function 
(Table 9) Instead of a Global Covariance Function (Table 1) 


Block 

Anom/m^g(mgal) Anom/m^g(mgal ) Anom/m^""g 

Current Procedures Current Procedures Dif f erences(mgal 
Global Cov. Function Rough Cov. Function Column 3- 
Terrestrial Data Set Terrestrial Data Set Column 2 

1105 

3. 4/5. 5 

3.2/10.0 

-0.2/4. 5 

1006 

13.0/5.3 

11.8/8.9 

-1.2/3. 6 

1007 

7. 9/3. 8 

8. 2/5.6 

0.3/1. 8 

1008 

8. 9/5. 3 

8. 8/9.0 

-0.1/3. 7 

1075 

21.8/3.8 

21.6/5.4 

-0.2/1. 6 

1076* 

3. 7/7. 5 

1.8/15.4 

-1.9/7. 9 

1077 

2. 4/4. 5 

3. 3/7. 5 

0. 9/3.0 

1078 

-4. 9/5. 6 

-3.9/10.2 

1.0/-4.4 

1142 

34.8/7.4 

28.7/14.4 

-6. 1/7.0 

1143 

14.1/12.3 

7.0/22.6 

-7.1/10.3 

1144* 

-1.3/4. 3 

0.0/6. 8 

1.3/2. 3 

1145 

-15.9/8.3 

-13.1/15.9 

2. 8/7. 4 

1208 

0.6/12.4 

0.1/22.2 

-0.5/9. 8 

1209* 

-20.3/7.4 

-21.3/12.3 

-1.0/4. 7 

1212 

-11.1/8.0 

-8.1/16.4 

3. 0/8. 4 

1211 

ft 

ft 

ft 

1270 

25.7/13.2 

13.9/23.1 

-11.8/9.9 

1271* 

-3. 2/7.6 

4.1/15.5 

7. 3/7. 9 

1272* 

-17.1/11.1 

-6.6/20.5 

10.5/9.4 

1273 

-16.0/8.7 

-10.3/18.- 

5. 7/9. 3 

1329 

15.3/5.4 

11.0/9.4 

-4. 3/4.0 

1330* 

-13.2/5.0 

-12.5/7.9 

0.7/2. 9 

1331 

ft 

ft 

ft 

1332 

-23.8/10.5 

-16.2/19.5 

7. 6/9.0 


ft No October 1979 terrestrial data exist for these blocks. 
* rough and very rough blocks 
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using current procedures with the global covariance function from 
Table 1 were copied from Table 6 , column 2 to Table 10, Column 
2. The differences between the two sets of predictions are given 
in Table 10, Column 4. 

It is seen that the anomaly differences given in Column 4 
are all within the standard errors of prediction given in Column 
2. In other words, there is no essential difference between the 
two sets of predicted anomalies. This result simply illustrates 
the conclusion already arrived at in Rapp (1977) that predicted 
anomalies arenot critically sensitive to the covariance function 
used. 


However, the significant effect of using different covariance 
functions is on the standard errors of predictions. Note the 
big differences in standard errors of predictions from the two 
sets, as given in Column 4. The much larger error of prediction 
associated with the use of a rough covariance function is caused 
by the enlarged variance and shortened correlation length for 
unchanged data spacing, implied by the use of the rough covariance 
function (see, for example, Sunkel (1981a)). Since a rough covar- 
iance function more truly represents the anomaly behavior over 
rough areas, it is to be expected that the standard errors of 
predictions given in Column 3 are more appropriate than those 
given in Column 2, in the case of rough blocks. For smooth blocks 
the standard errors of predictions should be those given in Column 
2 since as already concluded in Section 9 theses values are satis - 
factory. 

In conclusion, since the standard error of prediction is sensi- 
tive to the covariance functions used, therefore a more appropriate 
value for the standard error of prediction may be obtained by 
using a covariance function which truly represents the anomaly 
behavior in the area of interest. For example, in the case of 
current prediction procedures and even in the case of Method A 
discussed in Section 10, a global covariance function can be used 
for smooth and mild 300nm blocks which according to Table 4 comprise 
roughly 85% of the world, while a rough covariance function should 
be used for rough and very rough 300nm blocks which according 
to Table 4 comprise roughly 15% of the world. 

It is to be noted that the use of a covariance function that 
is appropriate for the anomaly field being considered will only 
solve the problem of the better estimation of the absolute error 
of prediction. Whether the absolute error of prediction will 
be large or small is another problem. To keep the error of predic- 
tion small, steps should be taken to smooth out a rough anomaly 
field in some way in order to obtain an anomaly field which has 
a smaller variance and a larger correlation length, (Sunkel 1981a, 
Figure 5). This will be the subject of succeeding sections. 
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12. Correlation Between 60nm anomalies and Elevations 


An examination of graphs of 60nm anomaly vs. elevation inside 
300nm blocks shows that there is a distinct linear trend between 
60nm anomalies and elevations. The linear trend is most pronounced 
in the case of 300nm blocks which are rough both in anomaly field 
and in terrain. This leads to the consideration of the possibility 
that 60nm elevations which happen to be fully known globally may 
contain enough information about 60nm anomalies that can be extracted 
to improve 60nm anomaly and eventually 300nm anomaly predictions. 

For the remainder of this report improvements to the 300nm 
anomaly prediction will be investigated with regards to the 
use of 60nm elevation information. In this section a study 
of the global behavior of 60nm anomalies and elevations and 
the global trend that exists between them will first be presented. 

The OSU combined altimeter /terrestrial data set of October 1979, being 
the most reliable and complete data set available at the time of 
the investigations, was used in the study to be presented. 


12.1 Global Behavior of 60nm Anomalies and Elevations 


Let X take on elevation values which are multiples of 
100 meters, and let Y take on anomaly values which are multiples 
of 5 mgals. Now define: 

Z xy (X,Y) = number of observed 60nm blocks in the world 
with the properties: 

X-50 < h < X+50 (meters) (56) 

Y-2.5 s Ag < Y+2.5 (mgals) , 
where P 

h : 60nm elevation 

Agp: 60nm anomaly. 

To give an example involving specific values, given X=-5500m 
and Y=-10 mgals, then Z is about 100, using the OSU combined 
altimeter /terrestrial data set of October 1979. A general picture 
of the plot of Zxy v s. X,Y is given in Figure 6, obtained from 
an analysis of the OSU combined altimeter /terrestrial data set 
of October 1979. 

Now let 

Z X (X) = number of observed 60nm blocks in the world with 
the property 

X-50 £ h < X+50 (meters), (57) 

without regard to Ag p . Zjj may be obtained form Z^y a s follows: 
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Figure 6: Histogram showing the Number of Observed 60nm Blocks 

- ^orld With 60nm Anomalies and Elevations 

Falling Within Specified 100-m Elevation and 
5-mgal Anomaly Intervals. 



(58) 


Z X (X) = I Z^CX.Y) 

Figure 7 shows a plot of vs. X , obtained from the Zxy 

shown in Figure 6. From Figure 7 it can be concluded that the 
60nm elevations of the world are predominantly clustered around 
two elevations, one at about -4200m and the other at about +300m. 

Again let 

Zy(Y) = number of observed 60nm blocks in the world with 
the properties: 

hp > -1200 (meters) (59) 

Y-2.5 < Agp < Y+2.5 (mgals), 


where -1200m is an arbitrarily selected X-value + that is found 
in between the two major peaks in Figure 6. Zy may be obtained 
from Zxy as follows: 

ZJ(Y) = l Z XY (X,Y> . ( 

X > -1200 ^ ' 

Figure 8 shows a plot of Zy vs. Y , obtained from the Zxy 
shown in Figure 6. Similarly, let 


Zy(Y) = number of observed 60nm blocks in the world with 
the properties: 


h„ < -1200 (meters 
Y-2.5 < gp < Y+2.5 (mgals). 


(61) 


Zy may be obtained from Zxy as follows: 

Zy(Y) = l Z xv (X,Y) 

X < -1200 


(62) 


Figure 9 shows a plot of Zy vs. Y , obtained from the Zxy 
shown in Figure 6. From Figure 8 and 9 it can be seen that 
as expected the 60nm anomalies of the world are clustered around 
zero, considering large samples (compare with Rapp (1977), Figure 
2 ). 


More detailed versions of Figures 8 and 9 may be obtained 
by considering 

Zy(X 0 ,Y) = number of observed 60nm blocks in the world 
with the properties: 

X 0 -50 < h < Xo+50 (meters 
Y-2.5 < Afep < Y+2.5 (mgals), 
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( 63 ) 



Figure 7 : Number of Observed 60nm Blocks in the World With 60nm Elevations 

_ Falling Within A Specified 100-m Elevation Interval. 


FREQUENCY COUNT 


- 160.00 - 96.00 - 32.00 32.00 

RNOMflL Y (MGRL) 


96.00 


160.00 


Number of Observed 60nm Blocks in the World with 
60nm Anomalies Falling With a Specified 5-mgal 
interval, Considering only 60nm Blocks with 60nm 
Elevations Greater than or Equal to -1200m. 



FREQUENCY COUNT (Z Y ) *10 

.00 76.00 152.00 228.00 304.00 380.00 



where X 0 is kept constant while Y is allowed to vary. Z Y 
may be obtained as a profile of Z XY (see Figure 6): 

Z y ( X o , Y ) = Z XY (X 0 , Y) . (64) 

Several different plots of Z Y vs. Y may be constructed given 
different Xo-values. It is found from actual plots that the 
plot Z Y vs. Y always tends to be bell-shaped as in Figures 
8 and 9. This fact will be the starting point of discussions 
in the next section. 


12.2 Global Trend Between 60nm Anomalies and Elevations 


12.2.1 Correlation Slope by Line Fitting 


O 

In the last section it was indicated that plots of Z Y 
vs. Y , which are profiles of Z XY (see Figure 6 and equation 
(64)), always tend to be bell-shaped. This means that 60nm 
anomalies Ag p occuring within a given 100-meter elevation 
interval 

X-50 ^ h p < X+50 (meters) 

tend to cluster around some average 60nm anomaly value. This 
average 60nm anomaly value can be expressed as (see Figure 6): 

Y Z XY (X,Y) • (Y) 

A'g(X) = - (mgal ) (65) 

l Z xy (X,Y) 

A plot of Ag vs. X is given in Figure 10, obtained from 
the Z XY shown in Figure 6. Figure 10 becomes a portrayal 
of the global trend that exists between 60nm anomalies and ele- 
vations. 

Figures 7 and 10 should be compared with each other. For 
convenience, the graph of the elevation counts shown in Figure 
7 was superimposed on the graph of average anomalies shown in 
Figure 10, with the resulting graph shown in Figure 11. In 
elevation ranges where there are large samples of cases occuring 
(see major peaks of elevation count graph) the average 60nm 
anomaly tends to zero (see horizontal intercepts of average 
anomaly graph). This corresponds to the tendency of positive 
and negative anomalies to cancel each other when a large sample 
of anomalies is considered. 

As smaller samples are considered a weak correlation between 
60nm anomalies and elevations tends to show up. Figure 11 shows 
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Figure 10 : Average of 60 nm Anomalies Occuring Withi 

A Specified 100-meter Elevation Interval 
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that as elevation goes away from the central value -4200m towards 
higher values up to about -1200m, anomaly tends to go away from 
the central value zero towards larger values up to about +15mgals 
on the average. As elevation gets lower and lower than -4200m 
the anomaly tends to get smaller and smaller than zero. This 
is just exactly the definition of positive correlation, which 
is the mutual tendency of two quantities which are random around 
their respective central values to be both larger or both smaller 
than their central values at the same time. It can be said 
that 60nm elevations which are less than about -1400m (see left 
peak in Figure 7) correlate positively with their corresponding 
60nm anomalies (see peak in Figure 9). This correlation is 

expressible as the slope of a line that fits the left branch 

of the graph in Figure 10 up to an elevation of about -1200m. 

Similarly, it can be said that the right peak in Figure 
7 correlates with the peak in Figure 8. Figure 11 shows that 
as the 60nm elevation becomes higher and higher than +300m, 
the 60nm anomaly tends to be larger and larger than zero. This 
corresponds to positive correlation between 60nm elevations 
greater than about +300m and their corresponding 60nm anomalies. 

This correlation is expressible as the slope of a line that 

fits the right branch of the graph in Figure 10 down to an elevation 

of about +300m. At elevations between roughly -1200m to +300m. 
Figure 10 exhibits a transition between the left and the right 
slopes of the total graph. 

Based on the slope of the graph in Figure 10 it seems 
sufficient to consider two correlation slopes, one for negative 
elevations and the other for positive elevations. A slope of 
+0.0024 mgal/m was found to fit the negative branch of the graph 
in Figure 10, and this slope then expresses the very weak positive 
correlation that globally exists between negative 60nm elevations 
and their corresponding 60nm anomalies. A slope of +0 . 0141mgal/m 
was found to fit the positive branch of the graph in Figure 
10, and this slope expresses a larger correlation between the 
positive 60nm elevations and their corresponding 60nm anomalies. 

In other words, in a global average sense, an increase of 1000 
meters in a negative 60nm elevation corresponds to an increase 
of about 2 mgals in the 60nm anomaly and an increase of 1000 
meters in a positive 60nm elevation corresponds to an increase 
of about 14 mgals in the 60nm anomaly. 


12.2.2 Correlation Slope by Ratio of Covariance Functions 


Correlation Slopes may also be computed, not by straight 
line fitting to graphs such as the one in Figure 10, but by 
taking the ratio of two covariance functions. 

Following section 5.1, numerical mean anomaly-elevation 
cross-covariances can be computed as follows: 
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cov ( Ag , Ah ) = 1 Ahk 

l A j A k 


( 66 ) 


Similarly, numerical mean elevation covariances can be computed 
as follows: 

covCAg.Ag) = IM A k ( 67 ) 

I A j A k 

Here, Ah is the mean elevation of the block, centered to the 
average of all mean elevations being considered. 

Now assume that 

Ag-j = b Ah . (68) 

J J 

where b is some constant that represents the correlation slope 
between Ag and Ah . Then, equation (66) may be written as: 

covC Ag , Ah ) ■ 1 A) A k b Ah . ) Ah k 

I A j A k 

= b - l. 4 i. A k ^U-A h Ji 

I A j A^. 

cov(Ag,Ah) = b cov(Ah,Ah) . (69) 


Equation (69) more readily follows from formulas for propagation 
of covariance functions through linear transformations such 
as equation (68) (see Moritz (1980), pp. 86-87). From equation 
(69) the correlation slope may be computed as the ratio: 


cov(Ag, Ah) 
cov (Ah , Ah) 


(70) 


In section 12.2.1 it was pointed out that negative 60nm 
elevations are correlated with their 60nm anomalies, with a 
global correlation slope of 0.0024 mgal/m. This correlation 
slope is the b-value that can also be computed from equation 
(70), if only those 60nm blocks with negative elevations are 
considered in the computation of the numerical covariance functions 
cov(Ag,Ah) and cov(Ah,Ah) . 

Using equations (66) and (67), along with conditions (14), 
(15), and (16), 60nm numerical covariance functions cov(Ag,Ah) 
and cov(Ah,Ah) were computed from the altimeter/terrestrial 
data set (Section 5.1). Only those 60nm blocks with negative 60nm 
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Table 11 


Covariance Functions cov(Ag,Ah) and cov(Ah,Ah) 
as Computed Based on Equations (66) and (67), 

Using only 60nm Blocks with Negative Elevation. 

Elevations are centered to their computed mean value -4120m 

Anomalies are centered to zero. 


deg. 

rad. 

cov( Ag, Ah) 
mgal-meter 

cov( Ah , Ah) 
meter 2 

0.0000 

0.0000 

8560 

3576553 

0.9803 

0.0171 

5782 

2041906 

1.6214 

0.0283 

4543 

1840109 

2.5079 

0.0438 

3623 

1556291 

3.4644 

0.0605 

3388 

1367771 

4.3780 

0.0764 

2953 

1217362 

5.3677 

0.0937 

2467 

1124613 

6.2205 

0.1086 

3434 

1117875 


Table 11A 

cov( Ag ,Ah) /cov( Ah ,Ah) 
(mgal/m) 


0.0033 

0.0029 

0.0025 

0.0023 

0.0025 

0.0024 

0.0022 

0.0031 


mean ratio: 0.00264 mgal/m 

std. dev. : 0.00039 mgal/m 
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elevations and known 60nm anomalies were considered in the computa- 
tions, with the 60nm elevations centered to their computed mean value 
value of -4120m. The results are shown in Table 11. 

In terms of the numerical functions, equation (70) means 
that b can be computed as the ratio of corresponding pairs 
of values from columns 3 and 4 of Table 11. The results of 
dividing the values in Table 11, Column 3 with the corresponding 
values from Table 11, Column 4 are shown in Table 11A. It is 
to be noted that the ratios are about the same, as they should 
be if Equation (68) is more or less correct. The mean ratio 
of +0 ,0026mgal/m agrees well with the correlation slope of +0.0024 
mgal/m obtained by fitting a straight line to the negative branch 
of the graph in Figure 10. Tables similar to Tables 11 and 
12 may be generated for the case of positive elevations and 
the resulting mean ratio compared with the value +0.0141 mgal/m 
obtained in Section 12.2.1 from line fitting. 


12.2.3 Anomaly Field Roughness vs. Terrain Roughness 


Since 60nm anomalies and elevations are correlated, though 
weakly in a global average sense, it would be expected that 
300nm blocks with a smoother anomaly field would tend to have 
a smoother terrain and those with rougher anomaly field would 
tend to have rougher terrain. 

Figure 12 shows the frequency count of 300nm blocks falling 
within the terrain roughness category of SMOOTH, MEDIUM ROUGH, 
ROUGH, and VERY ROUGH (see Table 4). Each bar is subdivided 
to show the distribution of 300nm blocks falling within the 
anomaly field roughness category of SMOOTH, MILD, ROUGH, and 
VERY ROUGH (Table 4). The bar chart indeed shows a weak global 
tendency for smoother anomaly field to be associated with smoother 
terrain and for rougher anomaly field to be associated with 
rougher terrain. 

Figure 13 shows the mean RMSCA for each terrain roughness 
category. Recall from Table 4 that RMSCA (RMS Centered Anomaly) 
is the root mean square of 60nm anomalies inside the 300nm block 
after the mean 60nm anomaly for the 300nm block has been subtracted 
from each 60nm anomaly. RMSCA being a measure of roughness of the 
local anomaly field inside the 300m block, shows that smoother 
anomaly field is associated with smoother terrain and rougher 
anomaly field is associated with rougher terrain. 


Figure 14 shows the mean RMSA for each terrain roughness 
category. Recall from Table 4 that RMSA (RMS Anomaly) is the 
root mean square of 60nm anomalies inside the 300nm block without 
removing the mean 60nm anomaly for she 300nm block. Figure 
14 shows the same conclusions as Figure 13, only this time the 
bars are about 5 mgals taller. 
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STATISTICAL ANALYSIS SYSTEM 
FREQUENCY BAR CHART 
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Figure 12 : Bar Chart Showing a) the Number of Observed 

300nm Blocks in the World with Terrain Roughness 
Falling into a Specified Category, and b) the 
Distribution of Anomaly Roughness Categories 
Occuring Within a Given Terrain Roughness Category. 
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STATISTICAL ANALYSIS SYSTEM 
BAR CHART OF MEANS 


RMSCA MEAN 


35 + 


30 + 


1 



***** 

1 



***** 

25 + 



***** 

1 



***** 

1 



***** 

1 



***** 

1 



***** 

20 + 



***** 

1 



***** 

1 



***** 

1 


***** 

***** 

1 


***** 

***** 

15 + 


***** 

***** 

1 


***** 

***** 

1 


***** 

***** 

1 


***** 

***** 

1 

' ***** 

***** 

***** 

10 + 

***** 

***** 

***** 

1 

***** 

***** 

***** 

1 

***** 

***** 

***** 

1 

***** 

***** 

** -.** 

1 

***** 

***** 

*:> .** 

5 + 

***** 

***** 

* * * * * 

1 

***** 

***** 

***** 

1 

***** 

***** 

***** 

1 

***** 

***** 

***** 

1 

***** 

***** 

***** 


**** 51 : 

***** 

***** 

***** 

***** 

***** 

***** 
***** 
***** 
***** 
***** 
***** 
***** 
***** 
* ** * 
***** 
***** 
***** 
***** 
***** 
***** 
***** 
***** 
***** 
***** 
***** 
***** 
***** 
***** 
***** 
***** 
***** 
***** 
***** 
***** 
***** 


SMTH TER MED RGH TER RGH TER VRY RCH TER 


ELEV VARTN 


Figure 13 : Average Value of RMS Centered 60nm Anomaly 

For a Specified Terrain Roughness Category. 
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STATISTICAL ANALYSIS 
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BAR CHART OF MEANS 


***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 




***** 

***** 



***** 

***** 



***** 

***** 



***** 

***** 



***** 

***** 



***** 

***** 



***** 

***** 



***** 

***** 



***** 

***** 


***** 

***** 

***** 


***** 

***** 

***** 


***** 

***** 

***** 


***** 

***** 

***** 


***** 

***** 

***** 


***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

if* rfC rfC 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 


***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

***** 

SMTH TER MED 

RCII TER 

RGH TER VRY 

RGH TER 


ELEV VARTN 


Average Value of RMS 60nm Anomaly 

For a Specified Terrain Roughness Category . 
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13. Removal of Local Trend Between 60nm Anomalies and Elevations 


It was shown in the last section that the correlation between 
60nm anomalies and elevations is weak in a global average sense. 
However, when local areas are considered, much larger correlations 
can be found. In particular, 300nm blocks with rough anomaly 
field, say RMSCA > 45mgals (see Table 4), exhibit much larger 
correlation slopes than the global values found in the last 
section. These correlation slopes are a form of local systematic 
effects on 60nm anomalies caused by local variations in 60nm 
elevations. It is therefore natural to investigate improvements 
of 300nm anomaly prediction through the use of the model of 
Least Squares Collocation with systematic parameters. 

Following equation (1), the vector of "observed” 60nm anomalies 
inside a 300nm block can be modeled as: 

Ag* = b Ah* + Ag* + n (71) 

where 

b : local correlation slope between 60nm anomalies and 

elevations inside the 300nm block 
Ah*: vector of 60nm elevations corresponding to Ag* , 

centered to the area-weighted average of the 25 60nm 
elevations inside the 300nm block; subset of the 
Ah in equation (74). 
bAh*: systematic part of Ag*. 

Ag*: signal part of Ag*. 

n: noise part of Ag*. 

The vector of 25 predicted 60nm anomalies would be (see equation 
( 6 )): 

Ag = b Ah + Ag g ^ (72) 

Various methods that can be used to estimate the parameter b 
will be described in Sections 14, 15, and 16. Assuming an esti- 
mate of b has been obtained, the signals Ag s can be solved 
using equations (4) and (5), which with changed notations become: 

Ag s - C Ag s Ag| <( fig|4g| + C „n>^ < Ag * - bAh *> • < 7 3> 

The vector of 25 60nm elevation differences can be computed 
from: 

Ah = h - h • I (74) 
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where 

h : vector of 25 60nm elevations 

h : area-weighted average of the 25 60nm elevations 

1 : column vector of 25 l's . 

The quantity E can be computed as: 

E = Lh ( 75 ) 

where L is the area-weighted averaging operator given in 
equation (34): 

L = L A“gAg * (76) 

Once the 25 60nm anomalies inside the 300nm block have been pre- 
dicted through equation (72), the 300nm anomaly can be computed 
as (see equations (33) and (34): 

Ag = L/fgAg Ag . 

Substituting equation (72) into the last equation: 

Ag = L AgAg (bAh +Ag s ) 

= L AgAg bAh + L AlAg Ag s * 

A'g = bLAh + L A - gAg Ag s (77) 

The first term on the right-hand side of equation (77) can be 
shown to equal zero: 

bLAh = bL(h -El) 

= bLh - bLEI 

= bE - bELl 

= bh - bh • 1 

= 0 . 

Therefore, equation (77) reduces to: 

Ag= L mg Ag s • (78) 

Equations (73) and (78) constitute the 300nm anomaly prediction 
equations under the model expressed by equation (71). 

As mentioned in section 11, in order to make the error of 
prediction smaller, steps have to be taken to smooth out a rough 
anomaly field in some way in order to obtain an anomaly field 
which has a smaller variance and a larger correlation length. 
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Using the model expressed by equation (71) is a step in this 
direction. In equation (73), the data vector is the smoothed 
(Ag* - bAh*) instead of the original vector Ag* used in the 
model Least Squares Collocation without systematic parameters 
(equation (54)) . The smaller the value of b , the less becomes 
the difference between the use of (Ag* - bAh*) compared with 
the use of Ag* as data vector in the prediction. As it will 
turn out in the next sections, the value of b for very rough 
blocks is large enough so that significant improvement in the 
prediction over rough blocks can be achieved by modeling b 
and using (Ag* - bAh*) instead of Ag* as data vector. 


14. Local Trend Removal by Least Squares Slope 


As shown in Figure 4, the test area over the Tonga Trench 
had good data coverage in the altimeter /terrestrial data set. 

Using this data set, Least Squares straight lines were fitted 
to plots of 60nm elevation vs. 60nm anomaly inside the 300nm 
test blocks. Equal weights were given to the known 60nm values. 
The slopes of the regression lines, which are actually the b- 
values, are given in Table 12. Also given are the characteristics 
of the 300nm blocks as to roughness of anomaly field as copied 
from Table 3. 

The rough 300nm blocks (marked with astrisk in Table 12) 
tend to exhibit much larger b-values than the smoother blocks. 
Since all the 60nm blocks of the test area have negative mean 
elevations, the b-values should be compared with the global 
average of 0.0024mgal/m found in Section 12.2. It is seen that 
the local b-values in Table 12 tend to be larger than the global 
average. 

Questions about the statistical significance of the determined 
values of b will be dealt with later, in Section 16. In 
the meantime, it is felt that the b-values in Table 12 are well 
determined because of the good coverage of the data set used, 
and it will be interesting at this point to go ahead and show 
prediction results using the computed b-values. 

Using equations (73) and (78), the altimeter /terrestrial 
subset described in Section 9, the global covariance function 
shown in Table 1, and the b-values from Table 12, predictions 
were done, with the results shown in Table 13. This method, 
which was Least Squares Collocation with the systematic parameter 
b obtained from Least Squares, will be called Method B. Only 
those 300nm test blocks which showed significant prediction 
differences are shown in Table 13. The other 300nm blocks had 
practically the same predicted results as compared with the 
predicted results from Method A (Section 10). Method A was 
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Table 12 



Slope b From Least Squares 
The Alt imeter /Terrestrial 

Using 

Set 


Slope b 


Block 

mgal /meter 

CHAR 

1005 

0.02375 

MLD 

1006 

0.01252 

MLD 

1007 

0.01464 

MLD 

1008 

0.01141 

MLD 

1075 

0.00616 

SMTH 

1076* 

0.03406 

VRY RGH 

1077 

0.03204 

MLD 

1078 

0.01432 

SMTH 

1142 

0.02444 

MLD 

1143 

0.01267 

MLD 

1144* 

0.03328 

VRY RGH 

1145 

-0.00689 

SMTH 

1208 

0.00194 

SMTH 

1209* 

0.02376 

VRY RGH 

1210 

-0.01507 

MLD 

1211 

0.01535 

SMTH 

1270 

0.00565 

SMTH 

1271* 

0.04299 

VRY RGH 

1272* 

0.04578 

RGH 

1273 

0.00027 

SMTH 

1329 

0.00429 

MLD 

1330* 

0.00163 

RGH 

1331 

-0.00971 

SMTH 

1332 

-0.00353 

SMTH 


♦rough and very rough blocks 
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Table 13 


Comparison of Predictions Without the Use of the 
b -Parameter (Method A) against Predictions With the 
use of the b-Parameter from Least Squares (Method B) . 

Altimeter /Terrestrial Subset Used. 


Block 

'true' anom 
alt /ter set 
(mgals) 

No. of 
known 
60nm 
blocks 
used 

anom/m^" 
(mgals) 
Method A 
alt /ter 
subset 

anom / m^-g 
(mgals; 
Method B 
alt /ter 
subset 

Anom 

Dif ferenc 
Method B- 
Method A 
(mgals) 

1076 

21.7 

13 

-4. 1/5. 5 

8. 7/5.5 

12.8 

1144 

10.3 

22 

6. 0/2.0 

9. 1/2.0 

3.1 

1209 

16.0 

14 

-11.9/3.7 

-0.2/3. 7 

11.7 

1271 

17.0 

12 

7. 3/5.7 

20.6/5.7 

13.3 

1272 

13.4 

5 

-2. 7/9. 7 

10.4/9.7 

13.1 


Table 14 

Comparison of Predictions Without the Use of the 
b-Parameter (Method A) against Predictions With the 
use of the b-parameter from Least Squares (Method B). 

Terrestrial Data Set Used. 



No. of 

anom /mA - g 

anom /mA-g 
(mgals) 

Anom 


known 

(mgals) 

Difference 


60nm 

Method A 

Method B 

Method B- 


blocks 

Oct 79 

Oct 79 

Method A 

Block 

used 

Terr Set 

Terr Set 

(mgals) 

1076 

13 

11.8/7.5 

5. 6/7. 5 

-6.2 

1144 

22 

1.6/4. 5 

6. 9/4. 5 

5.3 

1209 

14 

-6. 6/7. 8 

-2. 5/7.8 

4.1 

1271 

12 

0.4/7. 5 

17.8/7.5 

17.4 

1272 

5 

-13.5/11.4 

-3.8/11.4 

9.7 
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used for comparison since it is the most improved version of 
current prediction procedures at this point of the investigation. 

To be noted from Table 13 is the large improvement in the 
predicted values through the use of Method B. As compared with 
predicted values from Method A, the values from Method B are 
much closer to "true" values by the amounts given in the last 
column ("true" values come from Table 5). The values in the 
last column are all greater than the values of standard errors 
of predictions m*g , indicating that the difference between 
Method A and Method B is significant. There was no attempt 
to propagate the error of Least Squares estimation of b into 
the predicted values, hence the error of prediction for Method 
B was computed by the same formula used in Method A. 

Method B, together with the global covariance function in 
Table 1, was again applied to predict the 300nm test block 
anomalies, this time using the terrestrial data set. Results 
are given in Table 14. Again, it is indicated from Table 14 
that there could be significant differences between the use 
of Method A and Method B. 

The conclusion from this section is that it is promising 
to base the idea of improving 300nm anomaly predictions over 
rough areas on the use of the model Least Squares Collocation 
with the b-parameter. It is very encouraging that the blocks 
that had unreasonable predictions under current procedures (Table 
5) are the same blocks that had large improvements under Method 
B. The next section will give an idea of how much improvement 
could be attained through the use of a b-value that closely 
represents the true correlation between 60nm anomalies and eleva- 
tions inside the 300nm block. 


15. Local Trend Removal by Empirical Slope 


In the last section, a statistical determination of the b-value 
was performed inside the 300nm block. In actual production 
applications, there may not be enough known 60nm blocks inside 
the 300nm block from which to statistically determine the b- 
value. An alternative is to use a b-value that is known to 
reasonably represent the gradient of 60nm anomaly as a function 
of elevation in the area of interest. A possible source of 
such b-values would be any empirically derived regional averages 
of b-values. 

Woollard and Daugherty (1970) conducted extensive empirical 
studies of the magnitude of free air anomaly gradients that 
exist over various typical topographic forms in the oceans: 
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continental shelves, continental border trenches, island arcs and 
trenches, rises and ridges. The data analyzed were in the form 
of a voluminous amount of oceanic gravity traverse data across 
the various topographic forms. A complex situation was found 
that precluded the establishment of consolidated models that 
can define gravity gradients, given characteristic topographic 
forms. However, a large number of empirically determined grad- 
ients were given, covering a wide sample of the world's oceanic 
areas. A general conclusion was reached that oceanic free air 
anomaly gradients as a function of elevation might be expected 
to vary between zero and 0 . 069mgals/m, although frequent excep- 
tions can be encountered (Woolard and Daugherty (1970), p.218). 

An overall value that was found for the Tonga Trench area was 
0.06mgal/m (ibid., p.248). 

The studies of Woollard and Daugherty involved data spacings 
that are consistent with those of the present investigation, 
namely angular distances of about 1°. It is thus reasonable 
to use the value 0.06mgal/m as b-value in 300nm anomaly predic- 
tions over the Tonga Trench area. Prediction results, using 
this b-value in equations (73) and (78) and using the global 
covariance function from Table 1, are as shown in Table 15, 

Column 5. The data set used was the altimeter/terrestrial subset 
(Section 9). The method will be called Method C. Again, there was 
no propagation of an error in b into the predicted values, 
hence predicted accuracies are the same as in Method A and Method 
B (Table 13). 

For comparison, predictions were repeated using an arbitrary 
b-value of 0.04 mgal/m, with the results shown in Table 15, Column 
4. Alsogiven, in the last column of Table 15, are the prediction 
results and b-values from Method B of the last section. 

The empirically based b-value of 0.06 mgal/m yielded predic- 
tion results which are very close to "true" values. An exception 
to this is block 1271 for which a b-value of about 0.04 mgal/m 
seems to be the best one to use. It is seen from a comparison 
of predicted values associated with the different b-values of 
0.04 mgal/m, 0.06 mgal/m, and the Least Squares b-values, that 
the prediction is sensitive to the b-value used. This sensitivity 
of predictions to b-values indicates that it will be preferable 
to obtain the b-value using some uniform set of procedures that 
will yield a unique b-value; this will be the subject of next 
sections. 

The conclusion from this section is that, if no reliable 
statistical estimate of b-value can be made, then a b-value 
that is empirically known to apply over the topographic area 
of interest may be used for the removal of local trend between 
60nm anomalies and elevations. 
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Table J.5 


Comparison of Prediction Results from the Use of 
an Empirical b-value Against Prediction Results 
from the Use of the Least Squares b-value 
(Anom and m^-g in mgals; b in mgal/m). 


Block 

"True" 

Anom 

No. of 
known 
60nm 
Blocks 

anom/nu-% 
b=0.04 
alt /ter 
subset 

anom/niA- K 
b=0.06 
Method C 
alt /ter 
subset 

anom/niA-g 
Least Squares b 
Method B 
alt /ter 
subset 

1076 

21.7 

13 

12.2/5.5 

18.4/5.5 

8. 7/5. 5 , 

b=0. 03406 

1144 

10.3 

22 

10.0/2.0 

11.7/2.0 

9. 1/2.0 , 

b=0. 03328 

1209 

16.0 

14 

8. 3/3. 7 

17.6/3.7 

-0/2/3. 7 , 

b=0. 02376 

1271 

17.0 

12 

21.2/5.7 

25.9/5.7 

20.6/5.7 , 

b=0. 04299 

1272 

13.4 

5 

12.7/9.7 

14.5/9.7 

10.4/9.7 , 

b=0. 04578 


16. Local Trend Removal by Collocation Parameters 


16.1 Significance of Parameters 


It will be shown in this section that in the use of the model 
Least Squares Collocation with systematic parameters, the only 
parameter that should be modeled is the parameter b , and only 
when the parameter b satisfies some strict significance test 
should it be used in the actual prediction. 

First, consider the modeling of two systematic parameters. 

The vector of known 60nm anomalies, Ag* , is modeled as: 

Ag* = AX + Ag* + n (79) 

s 


where now, 

A = [I Ah*] and X T = [Ago b] (80) 


1 : column vector containing all l's. 

Ah* : vector of 60nm elevations corresponding to Ag* , centered 

to the area-weighted average of 60nm elevations inside 
the 300nm block. 

Ago : a constant shift of the 60nm anomaly values inside 

the 300nm block 

b : correlation slope between 60nm anomalies and elevations 

inside the 300nm block 
n : noise part of observations 
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Applying equation (2), the parameters Ago and b can be 
computed from: 


X = (A 


(C Ag* Ag* 


+ C ) 1 A) 


, -i 


nn 


(C 


' Ag *Ag* 


+ q, n ) 1 a g * (8i) 


The error covariance of the estimated parameters can be computed 
by applying equation (3): 


E XX ■ (A <C Ag*ig. + C nn rl A) ' 1 
s s 


(82) 


The variance of the predicted value of Ago is the element 
of E^x i n the (1,1) position: 


Ago 


E XX (1,1) ‘ 


(83) 


Similarly, the variance of the predicted value of b is: 


°b 2 E XX (2 ,2 ) * 


(84) 


Now, the t-test can be applied to measure the significance 
of the estimated parameters Ag 0 and b . This procedure 
involves the assumption of normally distributed residuals, 
but as Uotila (1967) pointed out the t-test is rather insensitive 
to moderate deviations from such an assumption. Moreover, 
it will be shown in the next section that the t-test indeed 
works as a device that can be used in a prediction program 
for the purpose of automatically deciding whether or not to 
use the parameter b in a certain 300nm anomaly prediction. 

The null hypotheses to be tested are: 

Hoi : Ago = 0 (85) 

H 02 : b =0 (86) 


It is to be recalled that the null hypothesis in statistical 
testing is a statement that is assumed to be true unless proven 
to be false beyond any "reasonable doubt." The amount of 
doubt that remains after the rejection of the null hypothesis 
is the type I error: the type I error is the probability 

that the rejected null hypothesis is in fact true. The user 
of a particular statistical test decides on the size of the 
type I error that can be tolerated. The size of the type 
I error is normally kept very small, since it is normally 
dangerous to reject a null hypothesis that is in fact true. 

In the case of the present report, it has already been 
shown that Method A, an improved version of current procedures, 
works satisfactorily over a great majority of 300nm blocks 
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of the world, and improvement is being sought only for rough 
blocks. Thus, in introducing improvements to Method A, it 
is necessary to disturb Method A as little as possible. Since 
Method A assumes Ago = 0 and b = 0 , therefore in order 
not to disturb Method A unnecessarily, the null hypotheses 
are chosen as they are in equations (85) and (86). 

The results of applying equation (81) to compute Ag 0 and 
b for the 300nm test blocks over the Tonga Trench Area (Section 
8) are shown in Table 16. The data set used was the altimeter/ter- 
restrial subset (obtained by considering only those altimeter/ter- 
restrial anomalies for which there was a corresponding known 
anomaly in the terrestrial data set; see Section 9), and 
the covariance function was the global covariance function 
given in Table 1. Also shown in Table 16 are the quantities 
I Ago / 0 | and |b/a b | , where a Ag 0 and a b were computed 
using equations (82), (83), and (847. Under the null hypothesis 
the quantities | Ag o / a Ag o I anc * |b/ofc,| are the quantities 
used in the t-test: 

t Ag 0 ( -’ m) = l Ago/a Ag 0 i (87) 

t b (a,m) = |b/a b | . (88) 

where 

a = probability 

m = degree of freedom 

(see, for example, Mikhail (1976), p.477). The degree of 
freedom m was taken to be 

m = no. of known 60nm blocks - 2 , (89) 

since there were two parameters being determined. 

Given t^ g0 (q,m) and the value of m , the probility a 
can be obtained from a t-table. The quantity a is the type 
I error if the null hypotehsis Hoi ' Ago = 0 ~ is rejected 
in favor of the predicted Ag 0 . Similarly, given the value 
of tt,(a,m) and the degree of freedom, the value of a obtained 
from a t-table is the type I error if the null hypothesis 
H o 2 : b=0 is rejected in favor of the predicted value of 

b . The type I errors associated with the Ago and b in 
columns 3 and 4 of Table 16 are shown in columns 7 and 8. 

As seen in column 7, the type I error associated with re- 
jecting the null hypothesis Ago=0 in favor of the 
predicted Ago in column 3 are all large. In other words, 
there is no reason to reject the null hypotehsis, namely 
Ago =0 . Indeed, there is no need to model a shift Ago : 
the approach being taken to improve predictions over rough 
areas is to smooth out the local anomaly field and the use 
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Table 16 


Testing for the Significance of b and A g 0 
in A g=A go +bAh+s+n 

(Least Squares Collocation Solution for b and Ago ; 
t-test to test Significance; 
Altimeter/Terrestrial subset used.) 


No .of 
Known 


60nm 

Blocks 

Block 

A gn 

(mgal ) 

b 

(mgal/m) 

| Ago /O ^g 0 1 

l b '<J 

t-test 

-Ag° 

t-test 
-b 

16 

1005 

8.3 

0.02070 

0.514 

2.314 

>25% 

1% 

16 

1006 

6.3 

0.00749 

0.394 

0.967 

>25 % 

20% 

22 

1007 

2.2 

0.02013 

0.143 

3.675 

>>25 % 

.05% 

16 

1008 

-2.5 

0.01985 

0.155 

1.985 

>>25% 

2.5 % 

23 

1075 

22.4 

0.00664 

1.416 

1.212 

10% 

10% 

13 

1076* 

8.5 

0.03586 

0.505 

8.019 

>25% 

<<.05% 

20 

1077 

3.2 

0.04676 

0.201 

2.225 

>>25% 

20% 

13 

1078 

-3.7 

0.01303 

0.221 

0.608 

>>25% 

25% 

11 

1142 

17.9 

0.03805 

1.051 

3.337 

15% 

.5% 

3 

1143 

21.4 

0.00675 

0.680 

0.216 

>>25% 

>>25% 

22 

1144* 

10.3 

0.04093 

0.654 

12.943 

25% 

<<.05% 

8 

1145 

8.0 

-0.00415 

0.461 

0.084 

> 25% 

>>25% 

3 

1208 

8.2 

0.00642 

0.236 

0.252 

>>25% 

>>25% • 

14 

1209* 

15.7 

0.03179 

0.970 

10.053 

20% 

<<.05% 

10 

1210 

11.3 

-0.04796 

0.659 

0.396 

25% 

>25% 

0 

1211 

ft 

ft 

ft 

ft 


ft 

2 

1270 

-3.5 

0.02000 

0.031 

0.396 

>>25% 

>>25% 

12 

1271* 

32.1 

0.05083 

1.835 

16.074 

5% 

<<.05% 

5 

1272 

41.2 

0.10485 

1.833 

6.054 

5% 

.5% 

7 

1273 

-0.1 

-0.00702 

0.005 

0.136 

>>25% 

>>25% 

19 

1329 

23.5 

0.00385 

1.458 

0.460 

10% 

>25% 

21 

1330 

2.0 

0.00703 

0.127 

1.572 

>>25% 

5% 

0 

1331 

ft 

if 

ft 

ft 

ft 

ft 

8 

1332 

5.7 

-0.03408 

0.314 

1.130 

>25% 

15% 
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of a Ago will not contribute to such smoothing; 
use of Ago burdens the solution for the parameter 
of the reduced degree of freedom. 


in fact the 
b because 


Column 8 of Table 16 shows that the type I errors associated 
with accepting the predicted b-values for very rough blocks 
(marked with asterisks) are very much smaller than the type T 
errors associated with accepting the predicted b-values for the 
other blocks. This suggests the possibility of using the t- 
test to automatically point to 300nm blocks where the parameter 
b needs to be modeled during predictions. For example, if the 
maximum type I error that can be tolerated is set to 0.05%, then 
only the b-values for very rough blocks will be found significant; 
for the smoother blocks, the b-value will continue to be taken 
as zero, since the null hypothesis H 0 2 : b=0 cannot be rejected 

for those blocks. 


To test the effect of not modeling Ag 0 , estimations were 
repeated, this time using: 

A = Ah* and X = b (90) 

in equations (79), (81), and (82). In this case, o\) 2 ~ E xx . 

The computed value of b and the t-test ratio | b / < 0 ^ | are 
shown in Table 17. Knowing t^(a ,m)= | b/Ab i and the degree of 
freedom m , the size of type I error a can be looked up from 
the t-table; ot-values are shown in the last column of Table 17. 

The b-values from Table 17 are very similar to the b-values 
from Table 16, except for those blocks with very few known 60nm 
blocks. The conclusions about which b-values are significant 
are unchanged in going from Table 16 to Table 17. 

In conclusions, only the parameter b should be modeled 
to improve predictions over rough areas; the shift Ago should 
continue to be modeled as Ago=0 . During the actual prediction, 
the estimated b-value should be used only when the type I error 
committed by using this b-value instead of b=0 is less than some 
very small tolerance. The next section will indicate the danger 
of using a b-value that is not highly significant, and will suggest 
a significance level for use in a production type of work. 


16.2 Recommended Significance Level 


For convenience, the equations for 300nm anomaly predictions 
using Least Squares Collocation with the systematic parameter 
b are summarized here. From equation (73), the observation 
vector is modeled as: 

Ag* = bAh* + Ag* + n . (91) 

o 
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Table 17 


Testing for the Significance of b in Ag = b^h+s+n; 

Least Squares Collocation to get b , t-test to test significance 

Altimeter /Terrestrial Subset from Section 9 used. 


No. of 


Known 


Block 

60nm 
Block s 

slope b 
mgal/m 

1 b/a b l 

t-test 

a 

1005 

16 

0.02049 

2.291 

2.5 % 

1006 

16 

0.00747 

0.964 

15% 

1007 

22 

0.02012 

3.673 

0.1% 

1008 

16 

0.01968 

1.968 

2.5% 

1075 

23 

0.00603 

1.101 

15% 

1076* 

13 

0.03557 

11.248 

<<0.05% 

1077 

20 

0.04673 

2.254 

2.5% 

1078 

13 

0.01294 

0.603 

25% 

1142 

11 

0.03995 

3.647 

0.4% 

1143 

3 

0.02255 

1.075 

25% 

1144* 

22 

0.04100 

12.965 

<<0.05% 

1145 

8 

-0.00400 

0. 081 

>>25% 

1208 

3 

0.01119 

0.722 

>25% 

1209* 

14 

0.03158 

9.986 

<<0.05% 

1210 

10 

-0.04143 

0.343 

>25% 

1211 

0 

if 

ft 

ft 

1270 

2 

0.01847 

1.847 

>>25% 

1271* 

12 

0.04938 

15.615 

<<0.05% 

1272 

5 

0.08960 

5.908 

0.5% 

1273 

7 

-0.00701 

0.136 

>>25% 

1329 

19 

0.00332 

0.397 

>25% 

1330 

21 

0.00704 

1.574 

7.5% 

1331 

0 

ft 

ft 

ft 

1332 

8 

-0.03219 

1.085 

20% 


#No October 1979 terrestrial data exist for these blocks 
♦blocks with very significant b-values; very rough blocks 
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Applying equation (2), the parameter b can be estimated from 
(see also equation (81)): 

b = <Ah» T (C AgtAg * + C nn ) -1 Ah*) -1 Ah* (C 4ggAg , + ^r'Ag* 

(92) 

Applying equation (3), the variance of the estimated parameter 
is (see also equation (82)): 

°b - < Ah * T < C Ag*Ag* + C nn>~' 4h * ) ' 1 ' (93) 

S s 

The signals at computation points can be solved using equation 
(73): 

Ag s = C Ag s Ag| (C Ag|Ag| + C nn ) 1 (Ag * " bAh+) . (94) 

The vector of 25 predicted 60nm anomalies is given by equation 
(72): 


Ag = bAh +Ag s . (95) 

Finally, the 300nm anomaly can be computed using equation (78): 

Sg = L AlTAg Ag s • (96) 

The variance of the predicted 300nm anomaly can be derived as 
follows: 

First, applying equation (7) and (8), the error covariance 
of the 25 predicted 60nm anomalies is: 

_ i p 

E AgAg = C ASsAg s " C Ag s Ag| (C Ag|Ag| + C nn ) C Ag s Ag| + 

+ (HAh* - Ah) a* (Ah* T hT - Ah T ) , (97) 

where 

H = C Ag s Ag* (C Ag|Ag| + C nn> 1 
Now apply equation (50): 

T 

E A”gA"g _ L A"gAg E AgAg L A _ gAg 

2 

The quantity E.g,-- is the same as m^g , the variance of the 
predicted 300nm anomaly. Substituting equation (97) into (99) 
and freely employing propagation of covariance matrices result 
in : 


(98) 

(99) 


- 64 - 


*"i T 

m A"g = C A"gA"g " C AgAg* (C Ag*Ag* + C nn ) C AgAg* + 
+ H Ah* ol Ah*T ST 

D 

where_ 

H = C A~gAg* (C Ag|Ag* + C nn )_1 * 


( 100 ) 

( 101 ) 


The linear transformation expressed by equation (96), together 
with the relations ^AgAg = ^ anc * LAh = 0 from section 13 
have been used in obtaining equations (100) and (101). Equation 
(100) allows for the propagation of the error of parameter estima- 
tion into the prediction results, through the presence of the 
third term on its right-hand side. 


A 300nm anomaly prediction procedure which will be called 
Method D is defined as follows: for each 300nm block, the fol- 

lowing steps are done: 


1. First compute b and cr b using equations (92) and (93). 

2. Form the ratio: 

t (a ,m) = — . 

a b 

3. Compute the degree of freedom: 

m = number of known 60nm blocks - 1 . 

4. Given m and a pre-selected value a (significance level, 
e.g., 0.005) look up the critical t-test variable to(a,m) 
from the t-table. 

5. If t(a,m) < to (a,m) set b = 0 and a b = 0 ; otherwise, 
ignore this step. 

6. Predict the 300nm anomaly using equations (94) and (96). 

7. Compute the variance of the predicted 300nm anomaly using 
equations (100) and (101). 


The presence of step 5 in Method D in effect causes Method D 
to reduce to Method A when the computed parameter b is not 
found to be significant at the selected significance level a . 

It is important to choose a proper value for a , in order to 
ensure that the already satisfactory performance of Method A 
over certain blocks will not be unnecessarily disturbed. The 
next two sets of predictions will point out a proper a to use. 

First, Method D was used to predict 300nm anomalies for 
the Tonga Trench area (Section 8) using a significance level 
of a = 0.005. The data set used was the altimeter/ terrestrial 
subset from Section 9, and the covariance function was the global 
function in Table 1. Prediction results are shown in Table 18. 
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Table 18 


Comparison of Prediction Results from Method D 
(i.e., Least Squares Colocation with the b-parameter 
Tested by t-test at 0.5% Significance Level) 
and from Method A (i.e., Least Squares Collocation 
Without Systematic Parameter) 



'True ' 

Anom/m r g 
(mgal ) 

Anom/m /"g 
(mgal 7 

Anom 



Anom(mgal ) 

Method A 

Method D 

Difference 

Slope 


Alt/Ter 

Alt/Ter 

Alt /Ter 

Method D- 

b 

Block 

Subset 

Subset 

Subset 

Method A 

(mgal /m) 

1005 

11.1 

8. 7/3. 5 

8. 7/3. 5 

0.0 

0.0 

1006 

9.2 

6. 5/3. 2 

6. 5/3. 2 

0.0 

0.0 

1007 

2.4 

0. 1/2.0 

0. 0/2.0 

-0.1 

0.00142 

1008 

-2.9 

-2. 3/3.1 

-2. 3/3.1 

0.0 

0.0 

1075 

19.4 

19.0/1.6 

19.0/1.6 

0.0 

0.0 

1076* 

21.7 

-4. 1/5. 5 

9.2/5. 7 

13.3 

0.03557 

1077 

5.3 

7.0/2. 5 

7. 0/2. 5 

0.0 

0.0 

1078 

-4.3 

-4. 0/3. 9 

-4. 0/3. 9 

0.0 

0.0 

1142 

23.4 

21.3/5.4 

12.6/6.0 

-8.7 

0.03995 

1143 

21.6 

13.3/11.0 

13.3/11.0 

0.0 

0.0 

1144* 

10.3 

6. 0/2.0 

9. 9/2.1 

3.9 

0.04100 

1145 

4.1 

6. 0/6.1 

6. 0/6.1 

0.0 

0.0 

1208 

16.5 

8.2/11.3 

8.2/11.3 

0.0 

0.0 

1209* 

16.0 

-11.9/3.7 

3. 6/4.0 

15.5 

0.3158 

1210 

17.9 

10.0/6.3 

10.0/6.3 

0.0 

0.0 

1211 

-0.5 

ft 

if 

ft 

ft 

1270 

22.2 

19.0/11.8 

19.0/11.8 

0.0 

0.0 

1271* 

17.0 

7.3/5. 7 

22.6/5.8 

15.3 

0.04938 

1272* 

13.4 

-2. 7/9. 7 

23.0/10.7 

25.7 

0.08960 

1273 

0.7 

-0. 2/8.0 

-0. 2/8.0 

0.0 

0.0 

1329 

27.1 

27.9/2.7 

27.9/2.7 

0.0 

0.0 

1330* 

2.4 

-4.3/2. 1 

-4.3/2. 1 

0.0 

0.0 

1331 

7.7 

ft 

If 

ft 

ft 

1332 

6.2 

-0. 2/8.1 

-0. 2/8.1 

0.0 

0.0 


#No October 1979 terrestrial data exists for these blocks 
♦Rough and very rough blocks. 
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For comparison, "true" values from Table 5 and the prediction 
results from Method A are also shown in Table 18. It is seen 
that over rough areas (marked with asterisk), Method D took effect 
yielding predicted values that are much closer to "true" values 
than the predicted values from Method A. Over the other blocks, 
it is desirable that Method D reduce to Method A since Method 
A already performs satisfactorily over those areas. Hlock 1142 
points up the danger of unnecessarily changing over from Method 
A to Method D: Method A already had a good prediction for block 

1142, and a switch over to Method D gave a worse prediction. 

To avoid an unnecessary switch over from Method A to Method 
D, a smaller value of the significance level was used: a = 0.0005 

The application of Method D was repeated, and the new results 
are shown in Table 19. This time, Method D took effect only 
on those blocks where Method A yields unreasonable predictions. 

For comparison of Method A and Method D (a = 0.0005), the errors 
of the two methods (predicted value minus true value from Table 
5) are the values given in Table 19. The errors should be viewed 
with proper regard to the corresponding standard errors of predic- 
tions. At this point of the investigations, the most improved 
prediction procedure, which is Method D, seems to be performing 
very well. 
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Table 19 


Comparison of Prediction Results from Method D 
(i.e., Least Squares Collocation with 
the b-parameter tested by t-test 
at 0.05% significance level) 
and from Method A 
(i.e., Least Squares Collocation 
without systematic parameter). 



Error An A -& 
(mgal r* 

Error/fa*-_ No 
(mgal) 

. of known 


60nm 



Method A 

Method D 

blocks 

, Slope b 

Block 

alt/ter subset 

alt/ter subset 

N 

(mgal /m) 

1005 

-2. 4/3. 5 

-2. 4/3. 5 

16 

0.0 

1006 

-2. 7/3. 2 

-2. 7/3. 2 

16 

0.0 

1007 

-2. 3/2.0 

-2. 3/2.0 

22 

0.0 

1008 

0.6/3. 1 

0.6/3. 1 

16 

0.0 

1075 

-0.4/1. 6 

-0.4/1. 6 

23 

0.0 

1076* 

-25.8/5.5 

-12.5/5.7 

13 

0.03557 

1077 

1.7/2. 5 

1.7/2. 5 

20 

0.0 

1078 

0.3/3. 9 

0. 3/3.9 

13 

0.0 

1142 

-2. 1/5. 4 

-2. 1/5. 4 

11 

0.0 

1143 

-8.3/11.0 

-8.3/11.0 

3 

0.0 

1144* 

-4. 3/2.0 

-0. 4/2.1 

22 

0.04100 

1145 

1. 9/6.1 

1. 9/6.1 

8 

0.0 

1208 

-8.3/11.3 

-8.3/11.3 

3 

0.0 

1209* 

-27.9/3.7 

-12.4/4.0 

14 

0.03158 

1210 

-7.9/6. 3 

-7. 9/6. 3 

10 

0.0 

1211 


ft 

0 

ft 

1270 

-3.2/11.8 

-3.2/11.8 

2 

0.0 

1271* 

-9. 7/5.7 

5. 6/5. 8 

12 

0.04938 

1272* 

-16.1/9.7 

-16.1/9.7 

5 

0.0 

1273 

-0. 9/8.0 

-0. 9/8.0 

7 

0.0 

1329 

0.8/2. 7 

0.8/2. 7 

19 

0.0 

1330* 

-6. 7/2.1 

-6. 7/2.1 

21 

0.0 

1331 

ff 

ff 

0 

if 

1332 

-6. 4/8.1 

-6. 4/8.1 

8 

0.0 

#No October 1979 terrestrial data exist for 

these 

blocks. 

♦rough 

and very rough blocks. 


- 
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17. Comparison of Methods of Trend Removal 


Table 20 is a compilation of important prediction results 
from Sections 14, 15 and 16. The blocks shown are typical 
blocks over which predictions with the parameter b should 
be applied. The different sets of prediction were all derived 
from the same altimeter /terrestrial subset defined in Section 9. 

For easier inter-comparison the values shown are errors of 
predictions ( error=predicted value - "true" value from Section 9), 
along with the standard errors of predictions, m^ . 

The values of m^-g from Method B and Method C are the 
same. These standard errors were computed using the same equa- 
tion used by Method A (Section 10), and consequently, there 
was no propagation made of any error in the b-value into the 
value of m^ - g . In contrast, the value of m^-g from Method D 
includes the propagated , the standard error of Least Squares 

Collocation estimation of b . It is seen that the effect 
of at, on mjfg is very small for the typical blocks shown. 

From an examination of the different sets of errors predic- 
tions, it is seen that the prediction is sensitive to the b- 
value used. It is therefore important to have a well-defined 
set of procedures for computing the b-value that will be used 
in the prediction. The recommended procedure here is Method D 
with a significance level of 0.05% To increase the reliability 
of b-values, it is proposed not to compute b-values for blocks 
with fewer than 4 known 60nm anomalies, the number 4 being 
a somewhat arbitrary number that we recommend based on experience 
with further numerical tests. In those blocks with at least 4 
known 60nm anomalies, the computed b-value should be tested for 
significance only when the b-value falls within reasonable limits; 
it is proposed to use |b| < 0.2mgal/m , since out of a total 
of 245 significant b-values that we have found in a global analysis 
of the OSU July 1981 terrestrial data set (Rapp, private com- 
munication) the b-values are mostly on the order of 0.04mgal/m , 
with very few b-values reaching 0.08 to 0.09mgal/m , and only 
one value reaching the maximum value O.lmgal/m . For those 
blocks where the parameter b will not be used in the prediction, 
the following should be set: b = 0 and = 0 , in order that 
the equations of Method D reduce to the equations of Method A. 

The main advantage of the procedures just described is that they 
can be used in a production type of work in order to reliably 
decide on which of Method A and Method B should be used for a 
particular situation. 

The b-values from Least Squares were obtained by straight 
line fitting to. conceptual plots of the 25 known 60nm altimeter/ 
terrestrial anomalies vs. elevations inside the 300nm blocks. 

Since the Least Squares method is a very common method 
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used for line fitting, the Least Squares b-value and the cor- 
responding predictions serve as reference values with which to 
compare the values from Method C and Method D. It is seen that 
Method D and Method B produce very similar b-values and prediction 
results. Method C and Method B have very different prediction 
results because of the difference in b-balues used; the b-value 
of 0.06mgal/m in Method C was empirically derived from gravity 
traverses across the Tonga Trench area, and have no direct connec- 
tion with the altimeter/terrestrial data values inside the individ- 
ual 300nm blocks. 

It is seen that the use of the empirical b-value 0.06mgal/m 
which applies over the Tonga Trench area produced very small errors 
of prediction. The predictions using the b-value of 0.04mgal/m 
were included to further illustrate the sensitivity of predictions 
to the b-value (compare the b-values and predictions from columns 
3 and 4, for example). An empirical b-value, if available 
for the characteristic topography of the 300nm block, can be 
used over rough 300nm blocks where no reasonable and statistically 
significant b-value can be found by Method D. The disadvantage 
of the empirical method is that the b-value is not unique, 
since there are no consolidated models that can uniquely define 
the gravity gradient given the character of the topography. 


18. Prediction Using Both Known Anomalies and Elevations 
As Observed Quantities 


In the previous sections, elevation data were used, considering 
them as error-free quantities. The b-parameter was computed 
under this assumption, and so 60nm elevations enter only in 
the design matrix A while the 60nm anomalies are the observed 
quantities (see Section 16). An alternative is to consider 
both 60nm anomalies and elevations as observed quantities, 
with their associated standard errors known from the data set 
being used. Anomaly predictions can then be done using standard 
Least Squares Collocation formulas. In this case, aside from the 
usual autocovariance matrix of 60nm anomalies, there will be a need 
for an autocovariance matrix of the signal part of 60nm elevations 
and also a crosscovariance matrix between the signal parts 
of 60nm anomalies and elevations. 


18.1 Prediction Equations 


Let 

Ag* : NX1 vector of observed 60nm anomalies inside 

the 300nm block 

Ah : 25X1 vector of observed 60nm elevations inside 

the 300nm block, centered to some mean elevation. 
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The observation vectors Ah and Ag* can be modeled as: 


Ag* = Ag* + n Ag 

Ah = Ah s + n^ h , 

which follow the form of equation (9), with: 


X = 

'Ag*' 

ll 

CO 

• * 

'Ag*' 

a 

ii 

" n Ag 


. Ah 


- A V 


- n A h . 


( 102 ) 

(103) 


(104) 


Equations (102) and (103) are in accordance with the model of 
Least Squares Collocation without systematic parameters. The 
vector of 25 predicted 60nm anomalies can be modeled as: 


Ag - Ag s , 

where Ag s is obtained by applying equation (10): 
Ag - Ag s « C sgsS , <C S , S , + C nn r> x . 
Substituting equations (104) into equation (106): 


(105) 


(106) 


Ag tC Ag„Ag| ;C AgsA h s ] 

o ^ 

25 x (N + 25) 


C Ag*Ag| + C n Ag n Ag ;C Ag|Ah s ’ 


'Ag* 

T 

^" C Ag^Ah ; C Ah s Ah s + Cn AhnAh 


.Ah . 


(Nf25) X (N-f-25) 


(N+25) X 1 
(107) 


The above equation, which predicts the 25 60nm anomalies inside 
the 300nm block using the models expressed by equations (102) 
and (103), will be referred to as Method E. 


The 300nm anomaly can now be computed using equations 
(33) and (34): 

= L AgAg AB • 


(108) 


The variance of the predicted Ag can be computed from the 
error covariance of Ag (see equation (49)): 


2 

Vg 


L AgAg E AgAg L AgAg 
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(109) 


Applying equation (11) to express E^ g .g > and employing propaga- 
tion of covariance matrices, equation (109) becomes: 



C A gA _ g 


^ C AgAg* C A"gAh s -* 


C Ag|Ag| + ° n A g n A g ’ C Ag|Ah s ' 

- i 

T “ 
C A“gA g| 

C T . q + C 

L AgsAh s ’ AhgAhg nAh n Ah_ 


T 

- C A-gAhs_ 


( 110 ) 


18.2 Covariance and Crosscovariance Functions and Matrices 


The computation of a numerical 60nm anomaly covariance function 
cov(Ag,Ag) is described in Section 5.1. The computation of 
a numerical 60nm elevation covariance function cov(Ah,Ah) and 
the computation of a numerical 60nm anomaly-elevation cross- 
covariance function cov(Ag,Ah) follow similar procedures as 
described in Section 5.1, and are explicitly given in Section 
12.2.2. An example of cov(Ag,Ag) is given in Table 1, while 
examples for cov(Ah,Ah) and cov(Ag,Ah) are given in Table 
11. Following exactly the same procedures described in Section 
5.2 for the generation of the basic anomaly covariance matrix 
CAgAg from cov(Ag,Ag), basic covariance matrices C A h A h and 
c AgAh ma y also be generated from cov(Ah,Ah) and cov(Ag,Ah). 

Then, following Section 5.3, other needed covariance matrices 
in equations (107) and (110) can. be derived by covariance propa- 
gation from the basic covariance matrices. For example, the 
C^-gAho appearing in equation (110) can be computed by area- 
weighted averaging columns of C^g^ • 


18.3 Prediction Results 


Equations (107), (108), and (110) were used to predict 
300nm anomalies for the 300nm test blocks. The method is called 
Method E. The first set of predictions are given in Table 21. 

The data set used was the altimeter /terrestrial subset. The 
covariance and crosscovariance functions used were those given 
in Tables 1 and 11. (The 60nm elevations in the test area are 
all negative, hence it is reasonable to use cov(Ag,Ah) and 
cov(Ah,Ah) from Table 11). In accordance with Table 11, the 
given 60nm anomalies were centered to the mean value for negative 
elevations: -4120m. 

Table 21 shows the difference between Method E and Method 
A. The values for Method A were copied from Table 8, column 
4. It is seen that Method E yields practically the same predicted 
anomalies and standard errors of prediction as Method A. This 
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Table 21 


Comparison of Prediction Results from the Use of Both Known Anomalies 
and Elevations as Observed Quantities (Method E, equation (102) and 
(105)) Against Prediction Results from the Use of only the Known 
Anomalies (Method A, Section 10). cov(Ag,Ag) from Table 1; 
cov(Ag,Ah) and cov(Ah,Ah) from Table 11; 
altimeter /terrestrial subset (Section 9) Used. 

(units: mgals). 


"True" Difference 


Block 

Anom 

Alt /Ter Set 

Anom/m^-g 
Method A 

Anom/m^-g 
Method E 

Method E- 
Method A 

1005 

11.1 

8. 7/3. 5 

8.8/3. 4 

0.1 

1006 

9.2 

6. 5/3.2 

6. 8/3.1 

0.3 

1007 

2.4 

0. 1/2.0 

0. 1/2.0 

0.0 

1008 

-2.9 

-2. 3/3.1 

-2. 3/3.1 

0.0 

1075 

19.4 

19.0/1.6 

19.0/1.6 

0.0 

1076* 

21.7 

-4. 1/5. 5 

-2. 1/5. 4 

2.0 

1077 

5.3 

7.0/2. 5 

6. 9/2. 5 

-0.1 

1078 

-4.3 

-4. 0/3. 9 

-4. 2/3. 8 

-0.2 

1142 

23.4 

21.3/5.4 

20.6/5.3 

-0.7 

1143 

21.6 

13.3/11.0 

13.8/10.8 

0.5 

1144* 

10.3 

6. 0/2.0 

6. 2/2.0 

0.2 

1145 

4.1 

6. 0/6.1 

5. 4/5. 9 

-0.6 

1208 

16.5 

8.2/11.3 

6.0/11.1 

-2.2 

1209* 

16.0 

-11.9/3.7 

-10.4/3.6 

1.5 

1210 

17.9 

10.0/6.3 

9. 5/6. 2 

-0.5 

1211 

-0. 5 

ft 

0.0/15.7 

# 

1270 

22.2 

19.0/11.8 

17.1/11.5 

-1.9 

1271* 

17.0 

7. 3/5. 7 

8. 6/5. 6 

1.3 

1272* 

13.4 

-2. 7/9. 7 

-3.2/9. 5 

-0.5 

1273 

0.7 

-0. 2/8.0 

-0.9/7. 8 

-0.7 

1329 

27.1 

27.9/2.7 

28.2/2.7 

0.3 

1330* 

2.4 

-4. 3/2.1 

-4. 2/2.1 

0.1 

1331 

7.7 

# 

-0.1/15.7 

# 

1332 

6.2 

-0.2/8. 1 

-1. 2/7.9 

-1.0 

#No October 1979 terrestrial anomaly data exist for these blocks. 

♦rough and very rough 

blocks. 
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is understandable because as shown in Table 11A, the functions 
cov(Ag,Ah) and cov(Ah,Ah) used in the predictions imply a very 
small anomaly-elevation correlation of 0.026mgal/m. 

A second set of predictions was obtained by repeating 
the first predictions, this time using the terrestrial data 
set. The results are shown in Table 22. For comparison, the 
results for Method A applied to the terrestrial data set 
are also shown. Again, there is essentially no difference between 
the results of Method A and Method E. 

To test the effect of using a drastically different set 
of covariance functions, cov(Ag,Ah) and cov(Ah,Ah) were first 
computed by equations (66) and (67) using Ag centered to zero 
(i.e., Ag used as they are) and Ah centered to the mean 
60nm elevation of individual 300nm blocks. The computed functions 
are shown in Table 23. Note the big difference between these 
functions and the corresponding functions in Table 11. 

A third set of predictions was done, differing from the 
first set only in that the covariance functions from Table 23 
were used instead of those from Table 11. The 60nm elevations 
were centered to the mean 60nm elevation inside the 300nm block, 
to obtain the Ah which was used as part of the data vector. 
Prediction results are shown in Table 24. Again, the results 
are to be compared with corresponding values from Method A. 

It is seen from Table 24 that, again, there is no substantial 
difference between Method A and Method E, even with the use 
of very different covariance functions from those used in Table 21. 

The conclusion from this section is that Method E, which 
uses both 60nm anomalies and elevations as observed quantities 
(see equation (102), does not seem to be able to extract a signif- 
icant amount of information about the 60nm anomaly from the 
known values of 60nm elevations. The effect is that equations 
(107) and (110) of Method E produce essentially the same results 
as equation (54) and (53) of Method A. 
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Table 22 


Comparison of Prediction Results from the Use of Both known Anomalies 
and Elevations as Observed Quantities (Method E, equations (107) 
and (110)) Against Prediction Results from the Use of Only the 
Known Anomalies (Method A, Section 10). cov(Ag,Ag) from Table 1; 
cov(Ag,Ah) and cov(Ah,Ah) from Table 11; terrestrial data set 
(Section 8) used. (units: mgals). 


Block 

anom/m^-g 
Method A 
Oct 79 
Terr Set 

Anom/m^g 
Method E 
Oct 79 
Terr Set 

Difference 
Method E- 
Method A 

1005 

4. 8/5. 5 

5. 3/5.4 

0.5 

1006 

12.3/5.3 

12.7/5.3 

0.4 

1007 

8. 3/3. 8 

8. 1/3. 8 

-0.2 

1008 

5. 5/5.3 

5. 3/5. 3 

-0.2 

1075 

18.6/3.7 

18.9/3.7 

0.3 

1076* 

11.8/7.5 

12.2/7.4 

0.4 

1077 

3. 7/4. 5 

3. 5/4. 5 

-0.2 

1078 

-4. 6/5. 6 

-5. 0/5. 6 

-0.4 

1142 

32.4/7.4 

31.8/7.3 

-0.6 

1143 

13.5/12.4 

14.4/12.2 

0.9 

1144* 

1.6/4. 5 

2.0/4. 5 

0.4 

1145 

-14.1/8.5 

-14.9/8.4 

-0.8 

1208 

0.5/12.8 

-0.9/12.6 

-1.4 

1209* 

-6. 6/7. 8 

-5. 6/7.7 

1.0 

1210 

-10.1/8.0 

-10.8/7.9 

-0.7 

1211 

# 

-3.1/15.7 

# 

1270 

23.5/13.4 

22.7/13.1 

-0.8 

1271* 

0.4/7. 5 

2. 2/7.4 

1.8 

1272* 

-13.5/11.4 

-14.8/11.2 

-1.3 

1273 

-15.4/9.0 

-16.4/8.8 

-1.0 

1329 

27.4/5.2 

29.1/5.0 

1.7 

1330* 

-15.1/4.8 

-15.2/4.7 

-0.1 

1331 


-3.1/15.7 

a 

1332 

-20.9/10.7 

-22.5/10.5 

-1.6 


#No October 1979 terrestrial anomaly data exist for these blocks 
♦rough and very rough blocks. 
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Table 23 


cov(Ag,Ah) and cov(Ah,Ah). h centered 
to 300nm block. Ag centered to zero. 


rp 

Deg. 

Rad. 

cov(Ah , Ah) 
meter**2 

No. of product 
pairs used 

0.0000 

0.0000 

535581 

41350 

0.9795 

0.0171 

346666 

55716 

1.6192 

0.0283 

174593 

104078 

2.5080 

0.0438 

-25945 

143208 

3.4656 

0.0605 

-185872 

113752 

4.3763 

0.0764 

-310360 

69580 

5.3610 

0.0936 

-421861 

9568 

5.2456 

0.1090 

-474506 

298 


Deg. 

Rad. 

cov(Ag, Ah) 
mgal-meter 

No. of product 
pairs used 

0.0000 

0.0000 

6714 

37814 

0.9795 

0.0171 

3596 

101875 

1.6204 

0.0283 

1499 

189887 

2.5084 

0.0438 

-596 

261425 

3.4662 

0.0605 

-1782 

208486 

4.3770 

0.0764 

-2424 

127758 

5.3625 

0.0936 

-3037 

17865 

6.2476 

0.1090 

-3842 

614 
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Table 24 


Comparison of Prediction Results from the Use of Both Known Anomalies 
and Elevations as Observed Quantities (Method E, equations (107) and 
(110)) Against Prediction Results from the Use of Only the Known 
Anomalies (Method A, Section 10). cov(Ag,Ag) from Table 1; 
cov(Ag,Ah) and cov(Ah,Ah) from Table 23; 
altimeter/terrestrial subset (Section 9) Used. 

(units: mgals). 


Block 

"True" Anom Anom/m^g 
alt/ter set Method A 

Anom/m^g 
Method E 

Difference 
Method E- 
Method A 

1005 

11.1 

8. 7/3. 5 

8. 6/3. 3 

-0.1 

1006 

9.2 

6. 5/3. 2 

6. 9/3.0 

0.4 

1007 

2.4 

0. 1/2.0 

0. 1/2.0 

0.0 

1008 

-2.9 

-2. 3/3.1 

-2. 1/3.0 

0.2 

1075 

19.4 

19.0/1.6 

19.0/1.6 

0.0 

1076* 

21.7 

-4. 1/5. 5 

-0. 5/5.1 

3.6 

1077 

5.3 

7.0/2. 5 

7. 1/2.4 

0.1 

1078 

-4.3 

-4. 0/3. 9 

-3.9/3. 7 

0.1 

1142 

23.4 

21.3/5.4 

19.8/5.0 

-1.5 

1143 

21.6 

13.3/11.0 

11.9/10.2 

-1.4 

1144* 

10.3 

6. 0/2.0 

6. 7/2.0 

0.7 

1145 

4.1 

6. 0/6.1 

6. 2/5. 7 

0.2 

1208 

16.5 

8.2/11.3 

2.5/10.5 

-5.7 

1209* 

16.0 

-11.9/3.7 

-8.6/3. 5 

3.3 

1210 

17.9 

10.0/6.3 

10.0/5.7 

0.0 

1211 

-0.5 

# 

0.0/16.0 

ft 

1270 

22.2 

19.0/11.8 

11.2/11.1 

-7.8 

1271* 

17.0 

7. 3/5. 7 

7. 7/5. 2 

0.4 

1272* 

13.4 

-2. 7/9. 7 

-2. 4/9.0 

0.3 

1273 

0.7 

-0. 2/8.0 

0.2/7. 4 

0.4 

1329 

27.1 

27.9/2.7 

28.0/2.7 

0.1 

1330* 

2.4 

-4. 3/2.1 

-4. 5/2.1 

-0.2 

1331 

7.7 

ft 

0.0/16.0 

If 

1332 

6.2 

-0. 2/8.1 

-1.1/7. 3 

-0.9 


#No October 1979 terrestrial anomaly data exist for these blocks, 
♦rough and very rough blocks. 
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19. Computer Implementation of Improved 300nm Anomaly 
Prediction Procedures 


As stated in the Introduction, the computer program supplied 
by Kaula (1966), which have had additions to it introduced by 
Rapp (1972), is the current operational program being used for 
the global prediction of 300nra anomalies from known values of 
l°x 1° mean free air anomalies inside the individual 300nm blocks. 

This program is called program F234. In order to implement the 
recommended improved procedures developed in this report, we have 
introduced modifications and additional coding to program F234. 

Call the improved program as program F234A. 

Program F234A retains the procedures of F234 regarding the 
computation of a global numerical covariance function of 60nm 
anomalies. This function is computed from a preliminary analysis 
of the entire global data set being used (Section 5.1). The formation 
of the basic covariance matrix C^g^g (Section 5.2) is also 
retained. The difference is in the manner of propagation of 
C AgAg to obtain ^A"gAg and C A - gA"g needed in the computation 
of tne standard error of prediction: whereas F234 uses straight 

averaging, program F234A uses area-weighted averaging from the 
elements of CAgAg (Section 5.3). As such, F234A computes a 
new set of CA - gA"g and C.-g^g for each 300nm block encountered 
in order to account for tne specific 60nm areas occuring in that 
300nm block. This area-weighted averaging is a consequence of 
the assumption that the 300nm anomaly should be computed as the 
area-weighted average of the 25 predicted 60 nm anomalies inside 
the 300nm block (see Section 5.3). 

The prediction procedures implemented in F234A are those 
of Method D as summarized in Section 16.2. The prediction technique 
used is therefore that of Least Squares Collocation with the system- 
atic parameter b , where b is the correlation slope between 
60nm anomalies and their corresponding elevations inside the 300nm 
block. The parameter b is analyzed for each 300nm block to 
be predicted. In order to control the conceivable occurrence 
of an unreasonable b-value that can pass the t-test (although 
we have not found this to be a problem even in a global prediction 
that we have conducted) it was decided to compute the b-value 
for the block only if there are at least 4 known 60nm anomalies 
inside the block, and moreover, to consider only b-values that 
satisfy |b| < 0.2mgal/m . A considered b-value is tested for 
significance using the t-test at 0.05% level of significance. 

The critical values to used for conducting the t-test at this 
significance level (Mikhail, 1976, p. 477) are input into the 
program through a DATA statement. The argument used to obtain 
the proper to - value is the degree of freedom, which is taken 
to be the number of known 60nm anomalies minus one. For a 300nm 
block in which no reasonable and significant b-value is found 
by the described procedures, the b-value and its standard error 
are set to zero and the computations continued. 
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Program F234A reads l°x 1° mean elevations as input data in 
addition to the usual l°x 1° mean anomalies and their standard 
errors. These l°x 1° elevations are formed into 60nm elevations 
just like the l°x 1° anomalies and their standard errors are formed 
into 60nm anomalies and propagated standard errors (see Section 
4 for the data structure used in the prediction). The 60nm eleva- 
tions enter the computations for b and , and the subsequent 

prediction equations. 

The original data vector Ag* is reduced to the mean elevation 
of the 300nm block, as in (Ag* - bAh*), for use in the prediction. 

All the 25 60nm anomalies of the 300nm block are predicted using 
the reduced data vector (see equation (94)). The 300nm anomaly 
is taken as the area-weighted average (equation (96)) of the 25 
predicted 60nm anomalies, using as weights the computed areas 
of the 60nm blocks. The predicted accuracy of the 300nm anomaly 
is computed using equations (100) and (101). 

As the 300nm blocks of the world are processed one after 
another, a record is kept of 300nm blocks which cannot be predicted 
because they do not have any l°x 1° anomaly data. These 300nm 
blocks are output at the end of all processing, in the format 
needed for direct input into the separate program of Rapp (Section 
7) which handles the prediction of 300nm blocks containing no 
l°x 1° anomaly data. It is proposed to retain Rapp's procedures 
for predicting 300nm blocks containing no l°x 1° anomaly data, 
namely, the prediction of the unknown 300nm anomaly from the 10 
closest 300nm anomalies which have been predicted beforehand (this 
time, these prior predictions can be done using the procedures 
of this report). The prediction from the 10 closest previously 
predicted 300nm anomalies involve much larger areas than a 300nm 
area that we have considered in this report, and therefore the 
use of the b-value (which applies to local areas like a 300nm 
block) has no direct applicability in this case. 
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20. Prediction of Two New Sets of 1654 300nm Anomalies 


20.1 Generation of the Two 1654 Sets 


In this section we will describe the different operational 
steps that were recently taken to predict two complete sets 
of 1654 300nm anomalies. The starting data set was the OSU 
terrestrial data set of December 1981, which contained 42,585 
l°x 1° mean free air anomalies (see Figure A.l of the Appendix) 
referred to the gravity formula of the Geodetic Reference System 
1967. Out of the 42,585 anomalies, there were 6,413 (see Figure 
A. 2 of the Appendix) which were coded as having been predicted 
by geophysical correlation techniques, and were therefore not 
based on actual gravity measurements. In accordance with Rapp 
(1972) these geophysical anomalies were given special consider- 
ation. That is, two sets of 1654 300nm anomalies were generated, 
one based on the original 42,585 l°x 1° anomaly data set, and 
the other based on the lesser 36,172 l°x 1° anomaly data set 
obtained after removing the 6,413 geophysically correlated 
anomalies from the original "42,585" data set. 

As a first step, the most recent OSU l°x 1° mean elevation 
data set was merged into the "42,585" data set. For the purposes 
of the new prediction program (program F234A of Section 19), 
l°x 1° elevations were merged into the data set even in the 
case of those l°x 1° squares which did not have the l°x 1° anomaly 
known (in this case, the l°x 1° anomaly was coded as 999 to 
denote it as being unknown). Since l°x 1° elevations are fully 
known globally, there were a total of 64,800 l°x 1° mean elevations 
that were merged into the fundamental data set. The use of 
l°x 1° mean elevation information is a feature of the improved 
procedures as implemented in program F234A. 

The "42,585" data set was actually a data base that contained 
much more information about each l°x 1° square than were needed 
in the predictions of this section; these extra information 
persisted in the merger data set of the last paragraph. The 
next step was therefore to pull out from the merger data set, 
only four pieces of needed information about each l°x 1° square 
of the world: 

1. the l°x 1° mean elevation 

2. the l°x 1° mean anomaly (if unknown, the anomaly was coded 
as 999) 

3. the standard deviation of the l°x 1° elevation 

4. the standard deviation of the l°x 1° anomaly (for an unknown 
anomaly, the standard deviation was coded as 99). 

In accordance with the housekeeping procedures of program F234A, 
each of the four sets of information given above was arranged 
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into 15°-latitude belts for input into F234A. Two separate 
input data sets (each arranged into 15°-latitude belts) were 
created for use in F234A. The first set (DEC81WGP) contained 
the complete 42,585 l°x 1° anomalies which included 6,413 geophys- 
ical anomalies, and the second set (DEC81NGP) contained the 
lesser 32,172 l°x 1° anomalies which excluded the 6,413 geophys- 
ical anomalies. 

Two separate sets of predictions using F234A were done, 
corresponding to the two input data sets. The two numerical 
60nm covariance functions that were generated at the start 
of each prediction set and were used in the predictions are 
given in Tables A.l and A. 2 of the Appendix. 

For DEC81WGP (with geophysical anomalies), 300nm anomaly 
predictions using F234A yielded predicted mean anomalies for 
1,542 300nm blocks. The remaining 112 300nm blocks which did 
not contain any l°x 1° anomaly data could not be predicted by 
F234A. However, the standard sequence numbers (see Section 
3), and the latitudes and longitudes of the centers of the 
112 unknown blocks, were output by F234A in a format consistent 
for input into the follow-up prediction programs which would 
be used to complete the prediction of the 1654 set. Predictions 
using DEC81NGP (no geophysical anomalies) yielded mean anomalies 
for 1,484 300nm blocks, and the sequence numbers and latitudes 
and longitudes of centers of the remaining 170 unknown 300nm 
blocks . 

The follow-up procedures used for completing each of the 
two sets of 1654 300nm anomalies were identical. Follow-up 
procedures were started by computing a numerical 300nm anomaly 
covariance function based on the partial set of predicted 300nm 
anomalies from program F234A. The method of generation of 
the covariances followed the global averaging procedures described 
in Section 5.1, this time using 300nm anomalies and areas instead 
of 60nm anomalies and areas. Using these global averaging 
procedures, covariances for ^-values from 0° to 180° in ip- 
steps of about 2° (=A^of Section 5.1) were generated. The 
covariances were then interpolated using the International 
Mathematical and Scientific Library (IMSL) routines IQHSCU 
and ICSEVU to form a regular table giving covariances for \p- 
values from 0° to 180° in uniform ^-steps of 0?5 (see Table 
A. 3 and Figure A. 3 of the Appendix). During actual predictions 
for the unknown 300nm anomalies, covariances for other if»- 
values not in the table were linearly interpolated from the 
values in the table. The ^-values were referred to the centers 
of the 300nm blocks. 

A program (programF342 of the OSU Department of Geodetic 
Science and Surveying) was used to predict the remaining unknown 
300nm anomalies using the following as input: 
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1. table giving 300nm anomaly covariances for ({/-values from 
0° to 180° in ({/-steps of 0?5 

2. the predicted 300nm anomalies and related information (lati- 
tude and longitude of block center and accuracy of predicted 
300nm anomaly) as output by program F234A. 

3. the sequence numbers, and the center latitudes and longitudes 
of the unknown 300nm blocks. 

The Unknown 300nm anomalies were individually predicted fr.om 
the 10 closest known 300nm anomalies, using Least Squares Colloca- 
tion with the numerical 300nm covariance function (Rapp, 1977). 
Prediction results were output in a format identical with 
the output of F234A, and as a final step these prediction results 
were merged with the predicted 300nm anomalies from F234A, forming 
a complete set of 1654 predicted 300nm anomalies. 


20.2 Comparison of 300nm Anomaly Sets 


Performing the above described follow-up procedures for 
both DEC81WGP (with geophysical anomalies) and DEC81NGP (no 
geophysical anomalies) predictions, the results were two complete 
sets of 1654 300nm anomalies , (the listing of the two sets are 
available upon request from The Ohio State University Department 
of Geodetic Science and Surveying). As part of operational 
procedures, comparisons were made involving these two 1654 sets 
and other previously existing 300nm anomaly sets. Table 25 
gives the description and selected statistics of the considered 
sets. 


Since the altimeter anomalies and the combination solution 
anomalies refer to the gravity formula of the Geodetic Reference 
System 1980, while the terrestrial anomalies refer to the gravity 
formula of the Geodetic Reference System 1967, the following 
transformation should be considered: 

Ag 1967 = Ag 1980 + 0,83 + 0,08 sin2< K m &a.;Ls) . 

In addition, an atmospheric mass correction <Sg^ should be 
considered, so that the final altimeter or combination solution 
anomaly that should be compared with the corresponding terrestrial 
anomaly would be: 

Ag' = Ag igg0 + 0.83 + 0.08 sin 2 <{> - 6g A (mgals) . 

In the last equation, the total correction to Ag^ggg turns 

out to be small (not more than 0.04mgal for an altimeter anomaly and 

often less than 0.5mgal for a combination solution anomaly) 

for all values of the latitude <{> , and this correction was 

neglected in subsequent comparisons. 
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Description and Statistics of Various 300nm Anomaly Sets 

(Statistics in mgals) 
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*(WGP means "with Geophysical", NGP "no geophysical" 


Table 26 shows some comparisons that may indicate the 
effect of including or not including l°x 1° geophysical anomalies 
in the generation of a 300nm anomaly set. There is an RMS 
difference of 3.9mgals between the set that includes geophysical 
anomalies (DEC81WGP)and the set that does not include geophys- 
ical anomalies (DEC81NGP), and this RMS difference is within 
the RMS standard deviation of the anomalies used in the compar- 
ison as given at the bottom of the table. A comparison of 
DEC81WGP and DEC81NGP with DEC81C0MB, also given in Table 26, 
shows that the set with geophysical anomalies has a smaller 
RMS difference with the combination solution than does the 
set with no geophysical anomalies (this is expected since the 
Combination Solution included geophysical anomalies as data — 
see Rapp, 1981). Overall, Table 26 indicates that DEC81WGP 
and DEC81NGP are consistent with each other and with DEC81 
COMB to within the RMS standard deviations of the anomalies 
compared, and that the DEC81WGP is nearer DEC81COMB than is 
DEC81NGP . 


Table 27 shows other comparisons of interest between 
various 300nm anomaly sets. The December 1981 and the October 
1979 terrestrial sets, which differ because of the incorporation 
of new l°x 1° data and the use of modified prediction procedures 
in the December 1981 set, have an RMS difference of 2.9mgals 
which is well within the RMS standard deviation of about 
6.5mgals (see Table 25) for each of the two sets. The other 
columns of Table 27 show that the December 1981 predictions 
compare slightly better with the November 1981 Seasat Altimeter 
set, than does the October 1979 predictions. Both the December 
1981 and the October 1979 sets are consistent with the altimeter 
set to within the RMS standard deviations of the anomaly 
sets compared. 
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Table 26 


Intercomparison of 300nm Anomaly Sets, Using Only those 300nm 
Anomalies such that the Standard Deviation is Less Than or 
Equal to lOmgals, and that the Corresponding Predicted 300nm 
Anomaly Values in DEC81WGP and DEC81NGP Actually Differ With 
Each Other. (See Table 25 for a Description of the 300nm 
Anomaly Sets.) 



DEC81WGP 

DEC81C0MB 

DEC81C0MB 


-DEC81NGP 

-DEC81WGP 

-DEC81NGP 


(mgals) 

(mgals) 

(mgals) 

Min. Difference 

-27.5 

-21.7 

-24 . 0 

Ave. Difference 

0.4 

-0.3 

0.1 

Max. Difference 

40.8 

27.0 

44.7 

RMS Difference 

3.9 

6.2 

7.2 

No. of Diff. >±3mgals 

1 

0 

1 

No. of Anom. Compared 

510 

510 

510 


RMS Std. Dev. of 300nm Anomalies used in: 


DEC81WGP - 4 . 7mgals 
DEC81NGP - 5 . 6mgals 
DEC81C0MB - 3.6mgals 


s 

g 

o 

o 

CO 


w 

o m 

•H rH 
h d 
d bo 
> £ 


CO 


<D 


X! 

E-i 


G 

(D 

<d m 
£ co 

p 

<D <D 
CQ H 
XI 
<D d 
O Eh 
G 

CD G 
U -H 


<D 

«H 


P <D 
•H Si 


•ri 

u 

<0 o 


Si 

P 


P 

0 


CQ 

<D 

Q 


cn 

p 
CQ CD 


o w 

•rH 


P >> 
m p 
•H ni 
P £ 
ri O 
P G 
W <3 









21. Summary and Conclusions 


In summarizing the studies conducted in this report, it is 
convenient to carry out the discussion in terms of three general 
types of anomaly field inside the 300nm block: 

Type I field: field which does not cause any roughness-related 

prediction problem; 

Type II field: field which causes prediction problems because 

of its roughness; 

Type III field: field obtained after removing from a type II 

field the linear effect of 60nm elevations on 
60nm anomalies. 


The main purpose of this report has been to improve the performance 
of current 300nm anomaly prediction procedures over a type II 
field, while at the same time not disturbing the satisfactory 
performance of current procedures over a type I field. 300nm 
blocks which have a rough anomaly field (see Table 4) and which 
contain relatively few observed 60nm blocks are the blocks with 
a type II field. 

Improvements which were not necessarily designed for type 
II fields were first introduced into the prediction procedures 
(Section 10). These were: (a) the prediction of all the 25 

60nm anomalies, unknown or known, inside the 300nm block; and 
(b) the use of area-weighted averaging as opposed to unweighted 
averaging of the 25 60nm anomalies to obtain the 300nm anomaly. 

The use of a covariance function that is representative of 
the variance and correlation length of the anomaly field of interest 
will result in a better estimation of the absolute accuracy of 
prediction. Whether this absolute accuracy of prediction will 
be high or low depends on the character of the field as to its 
variance and correlation length. To increase the absolute accuracy 
of prediction, a type II field must be smoothed in some way in 
order to obtain an anomaly field with lower variance and longer 
correlation length (see Sunkel (1981a)). Prediction should then 
be carried out from the smoothed field, then whatever systematic 
influences were removed prior to the prediction must be put back 
into the predicted value to obtain the final predicted value. 

With a view towards smoothing a type II field by removing 
the linear effect of 60nm anomalies on 60nm elevations, the cor- 
relation between 60nm anomalies and their elevations was first 
investigated in a global average sense. Two methods for computing 
an average correlation slope, one using line fitting to a graph 
of mean 60nm anomaly vs. 60nm elevation interval, and the other 
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using the ratio of two covariance function, were presented (Section 
12). Weak correlation slopes (b-values) of 0.0024mgal/m for negative 
60nm elevations and 0.0141mgal/m for positive 60nm elevations 
were found. Charts showing the global association between anomaly 
field roughness and terrain roughness were given in Section 12.2.3. 

In order to improve 300nm anomaly prediction over a type 
II field, it is proposed that the model of Least Squares Colloca- 
tion with systematic parameters be used, the parameter being the 
correlation slobe b between 60nm anomalies and their 60nm eleva- 
tions. Three possible sources of the value of the b-parameter 
were discussed. The b-value from Least Squares served as reference 
value, since Least Squares line fitting to data values is a very 
common method of computing regression slopes. The use of an existing 
empirical b-value, one that applies to the characteristic topo- 
graphy of the 300nm block, can be resorted to if there is an insuf- 
ficient number of known 60nm blocks inside the 300nm block to 
determine a reasonable and significant b-value. Finally, the 
b-value computed by Least Squares Collocation equation is the 
recommended b-value to use if it turns out to be reasonable and 
significant . 

To eliminate a priori unreliable b-values, it is proposed 
to compute the b-value only for blocks with at least 4 known 
60nm anomalies. Also, it is proposed to consider only b-values 
which have |b| < 0.2mgal/m as reasonable b-values. Of the b- 
values considered, only those b-values which are statistically 
significant under a t-test at 0.05% level of significance should 
be actually used in the prediction. These requirements appear 
satisfactory in practice, because the b-values that pass these 
requirements and are thus actually used in the prediction turn 
out to be the b-values of 300nm blocks with type II field. 

Some numerical results were also presented on the prediction 
of 300nm anomaly using both anomalies and elevations as observed 
quantities. It appears that no additional information about the 
300nm anomaly being predicted can be extracted using such predic- 
tion procedure. 

Remarks about other study areas connected with this report 
should be made at this point. This report considered the removal 
of local linear trend between 60nm anomalies and their 60nm eleva- 
tions, essentially producing a type III field from a type II field. 

A type III field is not necessarily a type I field, although a 
type III field is closer to a type I field than does a type II 
field. The reason for this is that there are factors other than 
the linear trend between 60nm anomalies and elevations that make 
a type II field different from.a type I field. Therefore, although 
the type III field is a significant step in the conversion from 
a type II field to a type I field, other data reduction steps 
can be studied. For example, Woollard and Daugherty (1970, p. 21) 
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state that there are geologic effects on the anomaly at sea 
on the order of 10 to 60mgals that are not related to topography. 

An overview of the effects on anomaly at sea caused by elevation 
and other geophysical phenomena may be found in Kearsley (1977), 
for example. 

As seen in Table 11, the improvement of 300nm anomaly 
prediction over a type II field could be a factor of two, 
through the use of the b-parameter. It will be of interest 
to study the feasibility of further improving the prediction 
through the use of the terrain correction. A recent report 
dealing with terrain correction computations has been published 
by Sunkel (1981b). 

The computer program developed for the operational imple- 
mentation of the recommended prediction procedures in this 
report is described in Section 19. Operational predictions 
were actually performed as described in Section 20, in order 
to generate two new sets of 1654 300nm anomalies, one involving 
l°x 1° geophysical anomalies and the other involving no geophys- 
ical anomalies. The two new sets incorporate the effect of 
new and updated anomalies in the starting l°x 1° data set, 
and the effect of modifications to prediction procedures as 
introduced in this report. Statistical comparisons show that 
the two sets (called the December 1981 terrestrial sets) are 
consistent with each other and with the November 1981 Altimeter 
set and the December 1981 Combination Solution set. The December 
1981 terrestrial set with geophysical anomalies agree better 
with the combination solution, than does the set without geophys- 
ical anomalies (this is expected since the combination solution 
included geophysical anomalies as data — see Rapp, 1981). Also 
the December 1981 terrestrial set with geophysical anomalies 
agree better with the altimeter derived set, than does the 
older October 1979 terrestrial data set with geophysical anomalies. 
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Appendix 


This appendix contains figures and tables referred to 
in Section 20, in relation to the prediction of complete 1654 
sets of 300nm anomalies. 
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Figure A.l Location of 42,585 l°xl° Anomalies of December 1981 Terrestrial Data Set 
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Location of 6,413 Geophysically Predicted l°xl° Anomalies 
Of December 1981 Terrestrial Data Set 






Table A.l 


60nm Covariance Function covCAg.Ag), Computed as Described 
in Section 5.1 Based on the December 1981 l°x 1° Terrestrial 
Data Set that Includes Geophysical Anomalies. Anomalies Refer 
to the Geodetic Reference System 1967. 


DEG 




RAD 


cov( Ag ,Ag)=C(i|0 
(mgal 2 ) 


No. of Product 
Pairs Used 


0.0000 

0.0000 

0.9798 

0.0171 

1.6232 

0.0283 

2.5085 

0.0438 

3.4657 

0.0605 

4.3780 

0.0764 

5.3729 

0.0938 

6.2380 

0.1089 


824.56 

30395 

478.39 

38394 

342.35 

68836 

251.40 

93364 

219.84 

73430 

205.80 

44933 

175.86 

5999 

155.43 

201 


Table A. 2 


60nm Covariance Function cov(Ag,Ag), Computed as Described in 
Section 5.1 Based on the December 1981 l°x 1° Terrestrial Data 
Set with Geophysical Anomalies Removed. Anomalies Refer to the 
Geodetic Reference System 1967. 


DEG 


RAD 

cov(Ag,Ag)=C(ifO 
(mgal 2 ) 

No. of Product 
Pairs Used 

0.0000 


0.0000 

860.80 

26600 

0.9802 


0.0171 

488.58 

32831 

1.6266 


0.0284 

336.09 

57724 

2.5091 


0.0438 

247.99 

77347 

3.4649 


0.0605 

221.59 

60515 

4.3796 


0.0764 

213.33 

36829 

5.3815 


0.0939 

162.94 

4792 

6.2427 


0.1090 

167.64 

177 
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300nm Anomaly Covariance Functions as Computed (see Section 20) 
Based on the "known" December 1981 300nm Anomalies. 
Anomalies Refer to the Geodetic Reference System 1967. 
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Figure A. 3 


Anomaly Covariance Function of the 
December 1981 300nm Anomalies. 

Anomalies Refer to the Geodetic Reference 
System 1967.. 

♦Geophysical Anomalies Included 
xGeophysical Anomalies Removed 



60.00 90.00 120.00 

SEPARATION (DEGREES) 


150.00 




